hep-th/ 030 1006 



Phases of Af = 1 Supersymmetric Gauge Theories 

and Matrices 



Freddy Cachazo, Nathan Seiberg and Edward Witten 
School of Natural Sciences, Institute for Advanced Study, Princeton NJ 08540 USA 



M = 1 supersymmetric U(N) gauge theory with adjoint matter $ and a polynomial super- 
potential Tr has been much studied recently. The classical theory has several vacua 
labeled by integers (Ni, N2, iV&), with the classical unbroken gauge group \~[iU(Ni). 
Quantum mechanically, each classical vacuum leads to Yli N{ different vacua. As the 
parameters of are varied, these vacua change in a continuous (and holomorphic) 
fashion. We find that vacua associated with (Nx, N2, N^) can be continuously trans- 
formed to vacua with (Ni,N2, ...,Nk), thus leading to a new kind of duality. Traditional 
order parameters, like the Wilson loop and 't Hooft loop, sometimes distinguish different 
phases. We also find phases that are not distinguished by conventional order parameters. 
The whole picture of the phase diagram is reminiscent of the phase diagram of M-theory. 
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1. Introduction 

Dynamics of four-dimensional supersymmetric gauge theories has proved to be a re- 
markably rich subject - one that hopefully will ultimately be important for the understand- 
ing of nature! Among other attractions of these theories, models with M = 1 supersym- 
metry exhibit many subtle properties, such as dynamical generation of mass, confinement 
of charge, and dynamical symmetry breaking, that are seen in the world of strong inter- 
actions and are still not fully understood (for a review, see e.g. |J). One example that 
has been much-studied in recent years is a U(N) gauge theory with a chiral superfield $ 
in the adjoint representation and a general single-trace superpotential Tr W(5>), for some 
polynomial function W that we will take to be of degree k + 1 and suitably generic. This 
model has been studied by its relation to geometric transitions and mirror symmetry in 
string theory and more recently by its surprising relation |||7]] to a bosonic matrix 
model - a zero- dimensional model with an iV x N matrix M and potential Tr W(M), 
where iV — ► oo. In this paper, we will investigate some additional fascinating properties of 
this example. 

Let us first explain the original question that motivated our investigation. In this 
model, if the function W'(x) has critical points a{, i = 1, . . . , k, then a classical vacuum 
is chosen by taking $ to be a diagonal matrix whose eigenvalues are the a,{, taken with 
multiplicity Nf, the N are any non- negative integers such that Yli^i = N. We let n be 
the number of choices of i with Ni > 0. Clearly, n < k. At the classical level, the choice of 
$ spontaneously breaks the gauge group from U (N) to Yli U (Ni) (where we include in the 
product only the positive Ni). Locally, this group is isomorphic to U(l) n x f] i SU(Ni). 
At low energies, the SU(Ni) gauge theories become confining and massive, so quantum 
mechanically, the gauge group that is actually observed at low energies is U(l) n . 

Since this is independent of Ni, the question now arises of whether it is possible 
quantum mechanically to distinguish in a precise way between vacua that have the same 
value of n but different iVj. To make this question clear, start with a choice of W(5>) such 
that the a, are far apart and the classical description of the $ field is a good approximation, 
and put the theory in a vacuum characterized by a given choice of integers N{. Now vary 
the parameters in possibly passing through a region in which are not far apart, 

but then return to a semiclassical region in which the are widely separated. In such a 
process, can one interpolate from a vacuum with one set of integers Ni to a vacuum with 
another set of integers iVj? 
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Let us provide a more systematic framework for this question. For each choice of 
W($), the theory under investigation has only a finite set of vacua. We will allow the 
parameters in W to vary, keeping only the leading coefficient (and the microscopic gauge 
coupling parameter of the underlying U(N) theory) fixed. The varying parameters describe 
a complex manifold T = C fc .0 As we vary the parameters in T, the vacua fit into a 
moduli space M. which because of supersymmetry is a complex manifold. M. is a union 
of components (irreducible complex sub manifolds) M. a which may possibly intersect at 
singularities. Each component M. a is a finite cover of T, and hence by going to infinity 
in T, we reach in each A4 a at least one semiclassical "end" with a definite set of iVj. Our 
question is now whether the same M. a can have different ends with different sets of iVj. 

To avoid undue suspense, let us assert that the answer to this question has turned 
out to be "yes." In considering the question, however, it soon becomes clear that one 
really should ask all the questions that are usually asked in string theory and M-theory. 
What are the components of M. and what distinguishes them? What semi-classical ends 
do they have and to what extent do they intersect at singularities? What is the physics at 
these singularities? (Some examples of singularities have been studied in ||.) When M. 
does have distinct components Ai a , what kind of natural "off-shell" interpolations can one 
make between them? Without going back to the underlying U(N) gauge theory (which 
really gives the most incisive definition of the whole problem), what kind of a low energy 
phenomenological Lagrangian or a "string field theory" can one find that describes all of 
the M. a at once? 

Smooth interpolations between classical limits with different gauge groups have been 
found in some other four-dimensional models. For example, in 0, in the context of M- 
theory compactification on a G2 manifold, an interpolation from SO (8 + 2n) to Sp(n) 
was argued, along with several other examples. This was generalized by Freedman ||10|| , 
who considered some additional examples, including models with SU(N) gauge symmetry 
broken by Wilson lines to Yl7=i SU(Ni) x C7(l) n_1 , and proposed interpolation between 
different sets of Ni in that case. Because the symmetry breaking pattern is the same, one 
may wonder if the model is actually equivalent by some string duality to the one considered 
in the present paper. One obstacle to finding such a duality is that it is not clear what the 



superpotential W($) would correspond to in the model considered in fpTOfj . At any rate, 



1 In this counting, we discard an irrelevant additive constant in W(^>) and set the highest 

9k „fc+i _ 1 , 
fe+i x fc+r 



power in W(x) to -8hx k+1 = T^rX k+1 



the present paper is devoted to studying a field theory problem by field theory methods; 
our interpolations and dualities are purely field theoretic and do not depend on stringy 
phenomena. 

It will become clear that the subject we are exploring is very rich, and in the present 
paper we will really only describe some general properties and explore some examples. 
Most likely, many interesting phenomena remain to be uncovered. 

Since it turns out that A4 does have many components, one basic question is how to 
distinguish them. One fundamental invariant that we have already identified is the rank n 
of the true low energy gauge group. In many physical theories, another important invariant 
comes from the realization of symmetries: which subgroup of the global symmetry group 
is unbroken in a given component of the moduli space of vacua? In the present paper, this 
criterion will be somewhat less useful, because for generic W there generally are no global 
symmetries; still, the realization of global symmetries, and more generally, transformations 
of $ that must also be taken to act on the parameters in T, will give some important 
information. Additional important information will come from a more subtle but still 
standard type of order parameter that is provided by confinement: for an external Wilson 
loop in a given representation of U(N), does one observe a Coulomb law or an area law? 

All these standard order parameters will play a major role in the present paper, but 
we will see that they do not suffice. We also will encounter a more subtle sort of order pa- 
rameter that only can be defined using supersymmetry. It may be characterized as follows. 
Suppose that a theory with J\f = 1 supersymmetry has r independent chiral superfields 
O s , s = 1, . . . , r.i By supersymmetry, the O s (or more precisely their expectation values) 
are holomorphic functions on A4. Now suppose that r exceeds the complex dimension of 
the moduli space Ai. Then as functions on any given component Ai a , the O s inevitably 
obey some algebraic equations; if these equations are different for different M. a , this gives 
an order parameter of sorts by which the M. a can be distinguished. We will find many 
examples of this situation, where the M. a cannot be distinguished by any conventional 
order parameter that we can see but can be distinguished by the chiral equations that are 
obeyed on them. When this occurs, we do not know if the a-dependence of these chiral 

2 When we say that they are independent, we mean that they obey no relations in the chiral 
ring of the theory that are independent of the parameters gk in the superpotential. Since we 
are working over the parameter space T, we are treating the gk as variables rather than as fixed 
complex numbers. 
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equations should be regarded as an explanation for the existence of different branches, or 
merely description of the phenomenon. 



Organization Of The Paper 

This paper is organized as follows. Section 2 is devoted to general considerations. In 
section 2.1, we describe what can be predicted about the moduli space of vacua from a 
"weak coupling" (actually weak gauge coupling) point of view in which the effects of the 
superpotential W($) are assumed to set in at much higher energies than the effects of the 
gauge coupling. We focus on understanding what can be said about the moduli space of 
vacua using the expected confinement of the low energy gauge theory. 



In section 2.2, following previous literature [f4|,|Tl||, we consider a strong (gauge) cou- 
pling point of view, in which the superpotential is considered as a small perturbation of 
a strongly coupled gauge theory with W = 0; this theory has Af = 2 supersymmetry and 
a vacuum structure determined in [|T2"| - |r4]| . In previous work, the strong coupling point of 
view was used to analyze the vacua with a given rank n of the unbroken symmetry group. 
That analysis assumed that the iVj defined above are all nonzero and that the degree of 
W is not too large. We generalize the analysis to remove these restrictions. (This analysis 
is completed in Appendix A.) 

In section 2.3, we review a construction originally described in 0] whereby, given 
vacua for the U(N) gauge theory with any given set of iVj, and any integer t > 1, one can 
construct vacua for the U(tN) gauge theory with N[ = tNj,. Comparing to our analysis 
of confinement in section 2.1, we will see that the vacua arising this way are precisely the 
confining vacua in which the confinement index (defined in section 2.1) is t. Conversely, all 
vacua with confinement index t arise by this construction. Since the theory with N- = tN^ 
also has vacua without confinement or with lower confinement index (as will be clear in 
section 2.1), there must be additional vacua with N- = tNf, they will be studied (in 
examples) in section 3. 

In section 2.4, we present some ways in which different branches M. a can arise and 
some constraints on the possible classical limits of a given branch. We also describe the 
simplest mechanism by which different branches may intersect. 

In section 3, we will explore the question of interpolating from one set of Ni to another. 
For this, we want to consider the simplest case in which n is as small as possible. For 
n = 1, the gauge group is completely unbroken at the classical level, so the question of 
interpolating between classical limits with different unbroken gauge groups does not arise. 
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So the first case is really that the rank of the low energy gauge group is n = 2 which 
(in the minimal case with all Ni non- vanishing) arises for a cubic superpotential, k = 2. 
Since we do not know how to make the analysis for general N, we will explore in detail 
the cases of N = 2, 3, 4, 5, and 6. For these values of iV and k = n = 2, we determine in 
detail all branches of the moduli space of vacua, recovering the confining vacua that were 
found in section 2.3 and describing the vacua without confinement. For the vacua without 
confinement, we do get the promised smooth continuations between different values of Ni, 
including interpolation from {N U N 2 ) = (1,3) to (2,2) for A = 4, from (1,4) to (2,3) 
for N = 5, and from (1,5) to (2,4) to (3,3) for N = 6. We conjecture that for any N, 
smooth interpolations are possible on the Coulomb branch from any (iVj, N 2 ) to any other 
(iVi,iV 2 ). 

In this detailed description for 2 < A < 6, we also see how branches with n = 2 meet 
branches with n = 1 at singularities. In addition to these intersections at finite points on 
the moduli space, we find for 2 < A < 6 that all branches meet the Coulomb branch at 
infinity, where the semiclassical description is valid. 

In section 4, we consider the opposite limit of large n. For n = N, that is all Ni = 1, 
there is a unique vacuum for each function VF($) and so our question about the different 
branches does not arise. So we consider the case that A — n is small but positive. Here, the 
simple strong coupling analysis of section 2.2 predicts that there must be at least N — n + 1 
branches of vacua. As far as we know, they are not distinguished by conventional order 
parameters. We show, as suggested above, that they can be distinguished by determining 
which holomorphic functions of the chiral superfields vanish on each branch. We do not 
know how best to interpret this result. 

In two appendices, we give more details of the strong gauge coupling analysis. In 
appendix A, we extend the proof that relates the matrix model curve to the M = 2 curve 
to W of high degree. In appendix B, we show how the generalized Konishi anomaly [fl5"l , |l6 



which was derived and used in [|17| arises in the strong gauge coupling regime where the 



elementary gauge fields are not visible. In appendix C we list some useful properties of 
Chebyshev polynomials. Appendix D defines a magnetic index v which characterizes the 
confinement in a theory with massless photons. It is related to the confinement index t of 
section 2.3 through A = tv. 



2. General Considerations 



2.1. Confinement 

Here, we will see what we can say about the U(N) gauge theory with adjoint superfield 
$ by using confinement as an order parameter. First we consider the matter directly in 
field theory, and then we re-examine the issues using the matrix model. 

We consider a classical vacuum in which Ni eigenvalues of $ are placed at the i th 
critical point of the superpotential W. To keep the exposition simple, we will assume 
that the Ni are all positive. We also assume that the critical points of W are distinct, 
so that all components of $ are massive classically. Supposing that the underlying U (N) 
gauge theory is weakly coupled at the scale set by those masses, the low energy physics is 
simply that of the pure supersymmetric gauge theory with classical gauge group Gf ow = 

wxn: =1 w). 

If Ni > 1, each of the SU(Ni) theories becomes confining at exponentially smaller 
energies. For any Ni > 1, the SU(Ni) theory has Ni vacua, each with a mass gap. The 
total number of vacua for fixed W and fixed Ni is hence Y\7=i ^i- If only W and Gf ow 
are specified, and not the individual number Ni of eigenvalues at the i th critical point, 
then the total number of vacua is larger, as one must allow for permutations of eigenvalues 
among the critical points. 

In the limit of small gauge coupling, the different SU (Ni) factors decouple and each 
of them confines. This confinement can be diagnosed by the possible area law of the 
Wilson loops Wi of the various SU(Ni). In the full theory, these Wi are not well-defined. 
Correspondingly, the full SU(N) theory is not necessarily confining even if the sub-theories 
are. (In this discussion, the U(l) factor in the underlying gauge group U(N) = U(l) x 
SU(N) is unimportant, as it decouples and does not contribute to the dynamics.) For a 
precise criterion for confinement, we have to use Wilson lines of the full SU(N) theory. 

To probe confinement, we consider a large Wilson loop in some representation 1Z of 
SU(N), and ask if it exhibits an area law or a Coulomb law. All that matters about 
TZ is how it transforms under the center of SU(N), which is Z^, generated by uj = 
diag(e 27rj/,Ar , e 2iri / N ^ . . . ; e 27 ™/^). The reason that only the center matters is that if TZ and 
TZ' are two irreducible representations that transform in the same way under the center, 
then an external charge in the representation TZ can combine with gluons to make an 
external charge in the representation TV , so one of these external charges is confined (has 
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infinite energy) if and only if the other is. Turning one representation into another by 
combining with gluons is a process that we may call electric screening. 

To give an example with any desired action of the center, we can simply take TZ to be 
the tensor product of r copies of the fundamental representation, for some r > 0. If W is 
a Wilson loop in the fundamental representation, the Wilson loop for the representation 
TZ is just W r . Only the value of r modulo N matters. 

If two representations TZ and TZ' are unconfined, then so is the tensor product TZ®TV . 
So the set of values of r for which there is no confinement is closed under addition. There 
is no area law for r = N (since r = iV is equivalent to r = 0) , so the smallest positive value 
of r for which there is no area law is always a divisor of N . We denote this number as t 
and call it the confinement index. If the theory is completely confining (all Wilson loops 
with non-trivial action of the center have an area law), then t = N, and if the theory is 
completely unconfining (no Wilson loop shows an area law), then t = 1. 

Let us now give a few examples of the behavior of the SU(N) theory. Suppose that 
SU(N) is broken classically to SU(N — 1) x £7(1). All vacua will be unconfining simply 
because there is no confinement in U(l). The fundamental representation of SU(N) has 
a component which is an SU(N — 1) singlet, and this component feels no confining gauge 
forces. When we evaluate the expectation value of a Wilson line in the fundamental 
representation TZ, 



the gauge field is effectively an (iV — 1) x (N — 1) matrix and a commuting lxl component 
associated with the U(l) factor (plus massive fields that will not influence the question 
of confinement). In computing the expectation value in (|2.1|) , assuming that C is very 
large compared to the scale at which SU(N — 1) becomes strongly coupled, the dominant 
contribution will come from the "bottom" component of TZ. This component does not 
interact with the SU(N — 1) gauge fields, so there is no area law. Since this component 
interacts with the massless U(l) gauge field, the Wilson loop ( |2.1| ) exhibits the behavior 
of a Coulomb phase rather than a Higgs phase. 

Computation Of Confinement Index In An Example 

For a simple example with a non-trivial confinement index, we will explore the case 
that SU(N) is broken to SU(Ni) x SU(N 2 ) x U(l). The fundamental representation of 
SU(N) decomposes as (N x , 1) © (1,N 2 ) under SU(Nx) x SU(N 2 ). (We ignore the U(l) 




(2.1) 
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charges, which will not affect confinement.) The tensor product of N\ copies of (Ni, 1) 
contains a singlet. So, as in the above example, the Wilson line W Nl will show no area 
law. Likewise, W N2 will show no area law. 

Let t' be the greatest common divisor of N\ and N2. Absence of confinement for 
r = N\ and for r = N2 implies that there is no confinement for r = l\N\ + I2N2 for any 
positive integers 1%, I2. We can pick l\ and I2 so that l\N\ + I2N2 is congruent mod N to 
t' . So the confinement index t can be no bigger than t' , and if electric screening, which we 
have considered so far, were the only mechanism, it would equal t' precisely. 

However, there is another mechanism, magnetic screening. Let us first recall 't Hooft's 
description of confined phases of 577 (N) gauge theory. One introduces an 't Hooft loop, 
which we will call H; H is constructed using a twist by an element of the center of SU(N). 
The general loop order parameter is W r H s , where r, s both take values from to N — 1 
(alternatively, their values only matter modulo N) . Massive phases are described by saying 
for which values of r, s there is no area law; one says that charges with this value of r and s 
are "condensed." (In some situations, for example a confining theory obtained by a small 
perturbation of an M = 2 theory, one can make the intuition behind this language precise 
In particular, there are N possible confining phases; in the r th such phase, W r H 
has no area law, or equivalently a charge with quantum numbers W r H is "condensed." 

M = 1 super Yang-Mills theory with gauge group SU(N) has iV confining vacua, with 
each type appearing precisely once. In fact, an adiabatic increase of the theta angle by 2n 
has the effect of increasing r by 1. 

So when we break the underlying SU (N) theory to 577 (iVi) x SU (N2) x U(l), we get 
a low energy theory that indeed has N1N2 vacua, but they are not equivalent. They can be 
distinguished by the type of confinement. For each pair n, T2 (with < < Ni — 1), there 
is one vacuum in which W^Hi and W 7 ^ 2 -^ are unconfined. (Here we write W\ and W2 
for Wilson lines in the fundamental representations of SU{N\) and SU(N2), respectively, 
and similarly Hi and H2 are the 't Hooft loops of the two factors of the low energy gauge 
group. As we remarked above, this makes sense only near the weak gauge coupling limit.) 

This theory also has 't Hooft-Polyakov magnetic monopoles. If N = 2, N\ = N2 = 1, 
we take the basic SU(2) monopole whose magnetic field at infinity is (in a unitary gauge) 
a multiple of 

For any N\, N2, and N, we just embed the 577(2) monopole in SU(N) in such a way 
that it does not fit entirely in either SU(Ni) or SU(N2) - we take the "1" in the upper 
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left to act in S77(ATi) and the in the lower right to act in SU(N 2 ). (We could also 

embed the monopole entirely in SU(Ni) or entirely in SU(N 2 ), but this would give nothing 
interesting for our present purposes.) 

Now suppose we probe the (ri, r 2 ) vacuum of the full SU (N) theory with the external 
Wilson loop W Tl ~ r2 . Although r\ — r 2 is not necessarily a multiple of t' , we claim that 
this Wilson loop will not show an area law. The external SU(N) charge v\ — r 2 can, after 
symmetry breaking to SU(Ni) x SU(N 2 ), be divided between the two factors in various 
ways, and confinement will be avoided if there is any way to divide the charge to make 
the energy finite. We simply put v\ units of charge in SU(Ni) and — r 2 units in SU(N 2 ). 
(The charges are only well-defined modulo Ni and N 2 .) So in other words, W ri ~ T2 can 
behave in the low energy theory as W = W[ 1 W 2 T2 . But the external Wilson line W can 
be screened by the spontaneous nucleation from the vacuum of the magnetic monopole 
described in the last paragraph. Since the monopole has magnetic charge 1 in SU(Ni) 
and —1 in SU(N 2 ), it can be represented as H\H 2 X . So in conjunction with a monopole, 
W can be represented in the low energy theory as H t W 2 r2 K 2 x , and has no area law 
since (by the definition of T\ and r 2 ), W^Hi has no area law in SU(Ni), and W^ r2 H^ 1 
(which is obtained from W2 2 H 2 by reversing the orientation of the loop) has no area law 
in SU(N 2 ). 

The net effect is that there is no area law for W Nl , for W N2 , or for W ri ~ T2 . Hence 
the confinement index t is at most the greatest common divisor of Ni, N 2 , and b = r± — r 2 . 
We claim that this is the correct value, since by now we have considered screening by all 
of the electric and magnetic objects that exist in the theory. 

It is easy to see why the confinement index depends not on r\ and r 2 separately but 
only on their difference b. In fact, under 9 — > 6 + 2n, one has Ti — > Ti + 1, for i = 1, 2, while 
external Wilson loops are unaffected. The confinement index must be invariant under this 
operation, so can only depend on the difference of the r^. 

General Case 

We can easily extend this analysis to the general case that U(N) is broken down to 
Fir=i U(Ni). Electric screening ensures that there is no area law for W N \ i = 1, . . . ,n. 
In addition, we must consider magnetic screening. The low energy theory has Yl7=i 
vacua, which are characterized by giving integers r^, i = 1, . . . ,n (with < < Ni — 1) 
such that in the low energy SU(Ni) theory, W^Hi has no area law. A magnetic monopole 
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can be shared between SU(Ni) and SU(N i+ i) in such a way that it has charges HiH^M 
It follows, just as in the case with n — 2 that was considered above, that the external 
Wilson line W ri ~ Ti+1 has no area law. (To make the argument, we simply place t% units of 
charge in SU(Ni) and — r i+ i units in SU(N i+ i) and then screen this configuration using 
the monopole.) 

So if we let 6j = — r^i, the confinement index t is now the greatest common divisor 
of the Ni and the hi. 

There are t different types of branch with confinement index t. For a branch of the 
u th type, the operators that do not have an area law are W f and W U H. 



Interpretation In The Matrix Model 

The matrix model is described Q by a complex curve S: 

y 2 m = W'(x) 2 + f n - 1 (x), 



(2.3) 



where / depends on the gauge coupling and the iVj. We assume for simplicity that k = n; 
i.e. no Ni vanishes. £ is a double cover of the x-plane. With W being of degree n, the right 
hand side of ( |2.3[ ) has 2n zeroes; the projection of £ to the complex x-plane is branched 
at these 2n points. 






Aj A 2 A 3 

Figure 1: Branch cuts, depicted by zigzags connecting the of , and cycles A% surrounding the 
cuts, are depicted here for n = 3. 



The situation can be described particularly simply for weak coupling, where / is small. 
The polynomial W has n roots Gtj, and for small /, the polynomial (W) 2 + f has, for each 
i, a pair of roots af that are near a^.B We connect the af with branch cuts, as in figure 
1, so that y m is a single- valued function away from these cuts. 



3 One can more generally share a monopole between SU(Ni) and SU(Nj) for any but 
this monopole leads to nothing new because its charges are a linear combination of the charges of 
monopoles considered in the text. 

4 When the parameters are such that the roots af are real, we ask that a\ > a~ . More 
generally, the symbols ± are just labels. 
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In [I7J, the operator- valued differential 



T(x) = Tr 



dx 



x - $ 



(2-4) 



was introduced. It was shown that over a cycle Ai surrounding the i th cut (also shown in 
figure 1), the period of T(x) was 



T. 



(2.5) 



Ai 



(To keep the formulas simple, in this paper we define the contour integral § to include a 
factor of l/2m.) In the present paper, we want to go away from the weak coupling limit, 
so we cannot assume that the 2n zeroes of (W') + f are neatly paired up. Likewise, we 
will not have a distinguished set of Ai cycles. So in developing the theory, we will have to 
include, along with the Ai, the other compact cycles with which the Ai will mix. There 
are n — 1 of these, sketched in figure 2. The new cycles, which we call Ci, roughly connect 
&t t° a i+i5 for z = 1, . . . , 7i — 1. We define 



T. 



(2.6) 



One might guess that since the Ni are integers, and will mix with the Cj under strong 
coupling monodromies, the q are also integers. We will prove this at the end of this 



section, using results from []17j • It turns out that the periods of T ( |2.6| ) are the integers Ci 
only on shell; i.e. only after all the equations of motion have been solved. Off shell, before 
solving the equations of motion of Si, the periods around the Ai cycles (|2.5|) are given by 
the integers Ni but the periods around the Cj cycles are not constrained. A more detailed 
discussion of this point will appear below and in a separate publication. 




Figure 2: Choice of compact cycles d for n = 3. The part of the contour depicted hy a dashed 
line is on the second sheet helow the cut. 
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Now let us recall the matrix model construction of the gauge theory effective super- 
potential. Here one uses another differential, which in the matrix model is 



t>{ \ 9m m dx 

R(x) = — Tr — : 

N x-M 



(2.7) 



where g m is the finite 't Hooft coupling of the matrix model. (|2.7|) has been interpreted in 



gauge theory in [17|. Its periods are 



R — Si 



A, 



1 dT 



[2.1 



i B . 2-ni dSi 

Here Bi are the noncompact cycles sketched in figure 3. The A\ and Bi are a canonical 
basis of the homology of the compact surface made by adding points at infinity to E; their 



intersection pairings are fl Aj = Bi fl Bj = 0, A^ D Bj 
compact cycles Ci introduced in the last paragraph by 



>i 3 - 



Ci — B i+ i — Bi 



The Bi are related to the 



(2.9) 



Furthermore, Si 



32 



^2-Tr WaiW" 1 is the gluino bilinear of U(Ni) (a more rigorous 



definition of Si is given in [0), and T is called the "prepotential." 





— • — 



Figure 3: Non-compact cycles Bi for n = 3. The bullets represent the two punctures of the 
Riemann surface (that is, the two points lying above x = oo). The cycle Bi passes through the 
i th cut and connects the two punctures, which are on different sheets. 
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To obtain the matrix model expression for the superpotential, one starts || with the 
fact that the model is dual to string theory on a certain non-compact Calabi-Yau threefold 
X. On the threefold, the superpotential is f x H A ft, where O is the covariantly 
constant holomorphic three- form of X, and the flux H has integral periods on compact 
cycles. By integrating out the extra dimensions, one reduces this expression to an integral 
over a Riemann surface (subsequently understood as the matrix model curve). In this 
reduction, H descends to a differential form To whose periods on compact cycles are 
integers, and Vt descends to the differential R. The superpotential becomes W = T A 

Evaluating this in the usual way in terms of contour integrals over one-cycles and 
Bi (and remembering our factor of l/2ni in the definition of a contour integral), we get 
the usual matrix model formula for the effective superpotential for the Sf. 

W ef f{Si) = 2my2(<f T <!> R- (f T I R) (2.10) 

• KJAi JB % JBi JAi J 

We will modify previous treatments only in not assuming that § B Tq is independent of i. 

According to the familiar analyses Q of the superpotential, § A Tq = A^. In addition, 
we set § c Tq = c^. Let us compare these to the periods of T as given in equations ( |2T5| ) 
and (|2.6|) . Since T and T have the same periods on the A-cycles, and the A-cycles will 
mix with the C-cycles under strong coupling monodromies, it is natural to conjecture that 
also their periods on the C-cycles are equal, Ci = c^. As will be explained in detail in a 
separate publication, this holds on-shell, that is after imposing the equations of motion of 
the Si. (Off-shell, the Ci periods of T depend on the Si, while T has integral periods even 
before imposing the Si equations of motion.) 

The integrals of Tq over the noncompact cycles Bi must then be as follows. § B ^ Tq = 
—tq for some complex number To, and more generally 

/ T = -TQ-bi, (2.11) 

JBi 

with 

i-l 

5 To, like H, is real; we have denoted it To because R is a differential of type (1,0), so only 
the (0, 1) part of To actually contributes to the integral. 
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So we get for the effective superpotential 



n n n 



'dSi 

i=l i=l i=2 



The derivation of (|2.13| ) does not require the knowledge that To = T on shell. 

The only addition that we are making to previous studies of this problem is to include 
the last term, proportional to the hi. Let us explain what this term does. The Veneziano- 



Yankielowicz superpotential for SU(M) supersymmetric gluodynamics is |18 



W(S) = S [log(A 3M /S M ) + M] . (2.14) 

The logarithm is only defined modulo 2ni, so W is not defined on the S'-plane but on 
an infinite cover thereof where for the time being we allow all possible branches of the 
logarithm. On this infinite cover, the number of critical points of W is M; the equation 
for a critical point is indeed 

log(A 3M /S M ) = 0, (2.15) 

and this equation implies that S = A 3M , an equation that has the familiar M roots. 
(Moreover, after picking S so S M = A 3M , we must, to obey ( |2.15| ), pick the correct branch 
of the logarithm, so each solution of S M = A 3M leads to only one vacuum.) Now consider, 



instead of ( |2.14|) , the superpotential 

W(S) = M[S log(A 3 /S) + S] . (2.16) 

Again, this superpotential is defined on an infinite cover of the S plane, but now on 
this infinite cover, the equation for a critical point is Mlog(A 3 /5') = 0, which implies 
log(A 3 /S) = 0. There is only one root, namely S = A 3 . 

What has happened? Mathematically, there are really M different ways to define an 
infinite cover of the S'-plane on which the Veneziano-Yankielowicz superpotential (|2.14|) is 
defined. The reason is that when the argument of S shifts by 2-7T, log(A 3 J S M ) shifts by 
—2itiM. We do not need to allow all branches of the logarithm to get a space on which W 
is defined; it suffices to consider 1/M of all of the branches. We could pick an arbitrary 
integer b and say that when A 3 /5' is positive, we want log(A 3M /5 ,M ) to have an imaginary 
part congruent to 2nib ( mod 2-KiM). So the Veneziano-Yankielowicz superpotential can 
be defined on each of M different spaces Sb, 6 = 0, . . . , M — 1, each of which is an infinite 
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cover of the S'-plane. Chiral symmetry permutes these, so physically we should include all 
of them, which is equivalent to including all branches of log(A 3M / S M ). 

The superpotential W(S) is equivalent to W(S) defined only on <Sq. It therefore only 
describes one of the vacua. To describe the others, in this language, we should explicitly 
include additional branches on which the superpotential is 

W b = M[S log(A 3 /S) +S]+ 2nibS. (2.17) 

The physical meaning of b is very simple. A superpotential term which is an imaginary 
multiple of S is just a shift in the theta angle. So including b just shifts the theta angle 
by 2rvb, which will rotate confinement (condensation of an 't Hooft loop) to oblique con- 
finement (condensation of a mixture of 't Hooft and Wilson loops) . Thus, for given 6, the 
vacuum that is realized has condensation of a 6-dependent combination of 't Hooft and 
Wilson loops W b H. 

Since in this paper we will be interested in studying U(M) gauge theories rather 
than SU (M) , we would like to extend (|2 . 1 7| ) to this case. As in |17| , we define S = 



S— (w a ) 2 . Here S is the trace over the SU (M) fields and it is taken to be an independent 
chiral superfield in the effective theory. w a is the field strength superfield of the U(l) C 
U(M). A fact which was useful in [17] is the decoupling of this U(l). This decoupling 
is implemented by considering the "superfield" S = S + if^ a w a — ij) 1 il) 2 M, where if) a is 
an auxiliary anticommuting Lorentz spinor. The first term in ( |2.17| ) is then written as 
jd 2 ^ [±S 2 log(A 3 AS) + fS 2 ] = M [siog(A 3 /<?) + s\ , which is independent of w a . (In 
fact, any expression of the form J d 2 ipH{S) is independent of w a .) The 2-KibS term in 
( ^.17| ) can also be written equally well in terms of S in (in U(M)) or S (in SU(M)), since 
b is an integer and abelian gauge theory is invariant under a 2-k shift in the theta angle. 

Now let us go back to the matrix model superpotential ( p2 ■ 13| ) , replacing M in the 
above discussion with any of the ATj. For small Si, dT/dSi ~ Si [log(Af /S'i) + l] . So our 
effective superpotential for small Si is (we drop the terms which involve two derivatives of 
the prepotential J 7 ) 

Weff(Si) = { 2 ™oSi + NiSi [log(A 3 /^) + l] + 2tu^) + O(SiSj). (2.18) 

i 

The first term is the bare coupling. The second term includes the one loop renormalization 
of the coupling constant and the strong SU(Ni) IR dynamics. The last term O(SiSj) 
represents the perturbative contribution of the high energy theory. Our new term 2nibiSi 
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arises also from the strong IR dynamics of SU(Ni). It should be independent of the photon 
of U(l) C U(Ni) and therefore the constants bi are quantized. 

From ( p.!8|) it is clear that the integers bi represent relative shifts by 2-rrbi of the theta 
angles of the various SU(Ni). More explicitly, relative to the theta angle in SU(Ni), the 
theta angle in SU (Nk) is greater by 2nbk- So if SU (iVi) has condensation of VF ri -ff,i then 
SU(Nk) has condensation of W Tk H, where — n = bk- (Our definitions are such that 
bi = 0.) 

Thus, the bk that we have defined in ( 2.12 ) are the same as the ones we introduced in 
our general discussion of confinement. So the confinement index t is the greatest common 
divisor of the Ni and the bk, or equivalently, the greatest common divisor of the N and 
Cj. Since N and q are a complete set of periods of T (integrated over compact cycles), we 
can describe this by saying that the confinement index is the greatest integer t such that 
all compact periods of T/t are integral. This gives, in the context of the matrix model, a 
manifestly "modular invariant" definition of the confinement index. 

In section 2.3, we will (following W) describe an operation that multiplies iV by an 
arbitrary positive integer t, also multiplying the individual Ni by t and multiplying T by 
t. This operation only generates confining vacua. In fact, it multiplies the confinement 
index by t. We will see that all confining vacua arise by applying this operation, starting 
with a non-confining vacuum with a smaller value of N. So for any given N, the "new" 
vacua that cannot be predicted based on what happened for smaller N are the "Coulomb" 
vacua, the ones without confinement. For every N, and every set of Ni, there are Coulomb 
vacua - for example, those vacua in which any of the bj is 1. 

Finally, we explain how to deduce from the results of |T7j that the c$ are integers. 
Instead of studying the present problem using the matrix model curve, we can use the 
Af = 2 curve 

y 2 = P N (xf - AK 2N . (2.19) 

(According to ||, these curves are closely related, with P^ — 4A 2N = Hj <r _ ri ((W') 2 + f); 
we will not need here to explicitly use this relation.) In appendix A of |17]], it was argued 
that T = (P' N /y)dx. This can alternatively be written 

T = —!—d(P N +y) = dlog(P N + y). (2.20) 

Pn + y 

Like any logarithmic derivative, T has integer periods (recall that we include a factor of 
2^ in the definition of <f). In other words, §T over a compact cycle is the change in 
log(Pjv + y) around that cycle, and this is an integer. 



The value of r\ depends on the bare coupling tq and on how H is defined. 
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2.2. Considerations based on strong gauge coupling 

The opposite of the weak gauge coupling approach of section 2.1 is a strong gauge 
coupling approach in which the superpotential is regarded as a small perturbation of the 
Af = 2 gauge dynamics that one would have if W = 0. This approach was developed in || 
by using methods of [|19] . Here we will review this analysis and extend it to allow some of 
the Ni to vanish and to allow superpotentials of any degree. 

We write the superpotential as 



k 



w = y^_ Tr ^ r +i_ r 22V) 

Z-^t y _1_ 1 V ' 



r=0 



r + 1 



In our present discussion, the superpotential is considered a small perturbation of an Af = 2 
U(N) theory. We interpret W as an effective superpotential on the Af = 2 Coulomb moduli 
space by replacing Tr $ r by its vev (Tr $ r ), regarded as a function on the moduli space. 
In order to look for vacua in which the low energy gauge group is U(l) n , we must extremize 
the superpotential ( |2.21| ) constrained to submanifolds of the Coulomb branch where N — n 
monopoles of the Af = 2 theory are massless (physically, the constraint is imposed because 
the massless monopoles get masses and vevs when W is turned on). The perturbation by 
W lifts all vacua of the Af = 2 Coulomb moduli space, except for a finite number which 
survive. 

It is convenient to introduce, for each Af = 2 vacuum, a classical N x N matrix $ c / 
such that (Tr $ r ) = Tr <E>^ for r = 1, . . . ,N. Moreover, we set u r = ^Tr $ r cl for all 
positive integers r. For r = 1,...,N, the u r are independent and are the usual order 
parameters of the Af = 2 theory. For r > N, both (Tr $ r ) and Tr $^ can be expressed 
in terms of the u r of r < N, but for r > 2iV, as shown in Appendix A of [ I7|] , they are 
unequal. 

The classical vacua of this theory are obtained by setting all eigenvalues of $ and $ c z 
equal to roots of W'(z) = ^2^. =0 g r z r '. In the present section, we will take the degree of 
the superpotential to be k + 1 < N, so that the u r that appear in the superpotential are 
independent and (Tr $ r ) = Tr In Appendix A, we generalize to arbitrary k. 

At a generic point on the Coulomb branch of the Af = 2 theory, the low energy gauge 
group is U(1) N . We want to study vacua in which the perturbation by W leaves only 
U(l) n gauge group at low energies, with n < k. This occurs if the remaining degrees of 
freedom become massive when W ^ because of the condensation of iV — n mutually local 
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magnetic monopoles. This can happen only at points where (at W = 0) the monopoles are 
massless. Monopole condensation for W ^ can be seen by including the N — n monopole 
hypermultiplets in the superpotential, 

N-n k 

W eff= Yl M i( u k)QiQi + Y,9rU r+1 . (2.22) 

1=1 r=0 

Here qi, qi are the monopole fields and Mi(uk) is the mass of I th monopole as a function 
of the Uk- In a supersymmetric vacuum, the variation of W e ff with respect to qi and qi 
is zero. However, for suitable g r , the variation with respect to the u^s does not allow 
qiqi to be zero. Therefore the mass of the monopoles has to vanish, i.e. M/((it&)) = for 
/ = l,...,N-n. 

The masses Mi are known to be equal to periods of a certain meromorphic one-form 
over some cycles of the M = 2 hyperelliptic curve, 

y 2 = P 2 {x) - AK 2N (2.23) 

where Pn(x) = det(x — $ c ;) is a polynomial of degree N. 

In 0] , it was shown that it is more convenient to use the fact that at points with N — n 
mutually local massless monopoles, the Af = 2 curve degenerates as follows: 

y 2 = P 2 N {x) - AK 2N = F 2n (x)H 2 N _ n (x). (2.24) 

This factorization is satisfied on an n-dimensional subspace of the Coulomb branch on 
which the superpotential should be extremized to find the points that preserve M = 1 
supersymmetry. 

The condition fl2.24j ) can easily be incorporated by means of Lagrange multipliers 0] . 
Using such a superpotential, it was shown in [|J that on shell and when the degree k of 
W'(x) is equal to n, the highest n + 1 coefficients of F 2n (a;) are given in terms of W'(x) 
as follows, 

F 2n (x) = \w'{x) 2 + Oix"- 1 ). (2.25) 

Assuming that W is given and the problem is to determine the M = 2 vacuum or equiv- 
alent^ Pn(x), (|2.25| ) is used as follows. (|2.25|) determines F 2n (x) in terms of n unknown 



coefficients. These are then determined, together with the desired Pjv, by asking for ex- 
istence of a polynomial Hn- u such that P^ — 4A 2N = F2 n H^_ n , as in ( |2.24|) . It is also 
often convenient to study the inverse problem: starting with an M = 2 vacuum with N — n 
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massless monopoles - in other words, a P/v that enables the factorization in ( P-24Q for 
suitable -ff/v_ n and Fin - one asks what superpotential W would lead to this vacuum. For 
k = n, this problem has a unique solution (modulo an irrelevant additive constant) since 
if F2n(x) is known, ( |2.25| ) determines W'(x). 

The main result of this subsection is to give the generalization of ( |2.25| ) to k > n, or 
in other words to the case in which some of the iVj are zero. We repeat the proof of ( 2.25Q 
given in ||, using a slightly different way to introduce the constraints. In Appendix A, we 
make a further generalization to arbitrary /c, dropping the restriction k < N — 1. 

Let us add to the superpotential constraints imposing the factorization (|2.24|) . It is 
useful to write Hn- u (x) = YliLi n ( x ~~ Pi)i i- e - Pi are the locations of the double roots of 
P^ - 4A 2N . Then we take 

w V* /r7r I Pn(x) -2e,A N f P N {x) - 2e,A N 
W -ff = 2^ 9ru r+1 + ^[L l j — — dx + B t j (x _^ )2 dx 



r=0 i=l 
k N—n 



r=Q i=l ^ J P% J \ Pi) / 



(2.26) 



Here Pn{x) is a function of u r . Li and B>i are Lagrange multipliers imposing the con- 
straints, 6i = ±1 and the contour of integration is a large curve around x = oo that 
encloses all p^s. We also have included a factor of -J^ in the definition of the symbol § in 
order to avoid cluttering of equations. 

The variation with respect to B>i leads to 

(The last step in ( |2.27[ ) follows from the fact that Pn(x) = det(x — <3> c ;) = ]^[^ (cc — 4> l c i), 
where 4> l cl are the eigenvalues of $ c z. Since Pn(Pi) ^ (given ( |2.24j) and the fact that 
Hjq- n {Pi) = 0), we learn that 

Tr 1 —— = (2.28) 

Pi - $cZ 

We can use either ( |2.27| ) or ( j2.28| ) to solve for pi in terms of the u r , but let us not do it 
now. Equation ( |2.27| ) or ( |2.28| ) will be used in Appendix B. 

Minimization of W e ff = W e f f(u r ,pi, Li, B>i; g r , A) leads to Wi ow = Wi ow (g r , A). How- 
ever, we will not attempt to carry this out. Instead, we will try to get information about 
Fm at the minimum by continuing to study the field equations. 
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Let us consider variations of (|2.26| ) with respect to pi 



'N- 



pn(x) ( t — — + 

i ^ (x -piY (x -PiY 



dx = 2B 



x - pi) 6 



dx = 



(2.29) 



where we used ( |2.27|) . Assuming that ( |2.24j) does not have any triple or higher roots, ( |2.29| ) 
implies that B>i = 0. 

Let us now consider variations with respect to u r . We use 



P N (x) = exp (IV \n(x - $ cZ )) = x N exp l-J^ ^ 



r=0 



(2.30) 



Since Pn(x) is a polynomial in x, ( |2.30| ) can be used to express u r with r > N in terms of 
u r with r < N by imposing the vanishing of negative powers of x. (See Appendix A.) We 
can also write (|2.30| ) as 

oo 

U r . 
X' 



Pn{x) 



x N exp( 



r=0 



(2.31) 



only in terms of the independent u r , those with r < iV.@ The derivative of ( f2.31| ) with 
respect to u r is 

dP N (x) _ \P N (x)] ^ 3 ^ 



8Ur 



Pn{x) 



X' 



We are now ready to differentiate W e ff with respect to u r : 



9r-l 



Pn(x) 



X' 



N-n 



— dx 

, x — Pi 



N- 



Mr J^ dx 

x r X — Pi 



(2.33) 



It is convenient to multiply ( |2.33j ) with z r 1 and to sum over r (z is inside the contour of 
integration) 



W'(z) = j 



N- 



x — z *■ — ' X — Pi 

i=l ^ % 



(2.34) 



We define the polynomial Q(x) in terms of 

La 



N- 



Q(x) 



y 

f-f x-pi H N - n (x) 
i=i 



(2.35) 



7 The symbol [A]+ denotes the polynomial part of a Laurent series A. 
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dx. (2.36) 



and write ( |2.34|) as 

W'(z) = t U Q{X ) . Pn[x) 
J H N - n (x) x - z 

Since W'(z) is a polynomial of degree k, the polynomial Q(x) is of degree k — n, and 
we will denote it as Qk- n (its definition (|2.35|) determined it as a polynomial of degree 
N-n-l>k-n). 

Finally, notice from (|2.24|) that, 



Pn(x) = v / F 2n (x)H N _ n (x) + 0{ 



x N ) 



So we can freely replace Pn(x) by a/ F 2n (x)HN- n (x) in ( |2.36|) to get 



W\z)=l Qk - n{x) ^ Jx) dx. (2.37) 
/ x-z 

We can easily recognize ( |2.37|) if it is written 

W '( z ) = I ^Mdx (2.38) 
J x-z 

with 

y 2 m (x)=F 2n (x)Ql_ n (x). (2.39) 

( p.39|) is the matrix model curve (for a case in which only n of a possible k cuts have 
"opened up" ) , and ( |2.38|) is the corresponding "equation of motion" for the matrix model 
eigenvalues. 

From ( |2.38j ) we conclude that (|2.25|) has to be generalized to 



F 2n (x)Ql_ n (x) = W' k {xf + 0{x k ~ l ). (2.40) 



In the case when k = n, we recover (|2.25|) with Q = g n . ( |2.40| ) can be used just like ( |2.25| ) 



to determine the M = 2 vacuum if W is given. With W known, ( |2.40| ) determines F 2n -, 
assuming its leading coefficient is normalized to 1, in terms of n parameters, which are 
then determined, along with Pjv and H^- n , by requiring the factorization ( |2.24j) . For the 
inverse problem, things are a little different. If we start with an JV = 2 vacuum, so that 
Pni F 2n , and Hn- u are known, and we want to constrain W ', then ( |2.40| ) actually leaves 



k — n free parameters in the determination of W. In the classical limit, the meaning of 
this is that if all eigenvalues of $ are placed at n zeroes of W, then the expectation values 
(Tr $ r ) are entirely independent of the other n — k zeroes. 
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In appendix A, we prove that the generalization of ( |2.40| ) for arbitrarily large k is 



F 2n (x)Ql_ n (x) = W' k {xf + 0{x k ~ x ) (2.41) 

where 

Qk-n(x) = V k -N{x)H N - n {x) + Q N - n - 1 (x) 

and Vk-N{x) is a polynomial of degree k — N whose coefficients are new Lagrange multi- 
pliers. 

Note that (|2~40|) and fl2~4lD have the same structure. Moreover, they provide enough 



equations, as in the case of ( |2.25|) , to determine the vacuum if W is given, or to determine 
W modulo a physically sensible redundancy in the case of the inverse problem. 

It is interesting to note that this result is consistent with the fact that the matrix 
model does not have information about the rank of the field theory gauge group U(N). It 
only knows about the number of cuts, i.e. the number of non-zero N{. 

See Appendix A for the generalization to superpotentials of arbitrary degree. 

2.3. The Multiplication Map And The Confinement Index 

In |], a construction was described that, for any positive integer t, maps vacua of 
the U(N) theory with a given superpotential W to vacua of the U(tN) theory with the 
same superpotential. This operation maps classical limits with unbroken Yli=i to 
classical limits with unbroken JX =1 U(tNi) and so in particular leaves fixed the rank n of 
the low energy gauge group. We now review this construction and show that it multiplies 
the confinement index by t and that all vacua of U(tN) with confinement index t arise 
from the Coulomb vacua of U(N) under this multiplication map. This implies in particular 
that, for any N, the only "new" vacua, which cannot be deduced from a U(N/t) vacuum 
(for some t) via the multiplication map, are the Coulomb vacua. 

A key role is played in the analysis by the polynomials |2(]] that describe those M = 2 



vacua of the U(t) theory that have the maximum number of massless monopoles, namely 
t — 1. These are the vacua that survive when the M = 2 theory is perturbed by a quadratic 
superpotential W = |Tr $ 2 ; this perturbation removes $ from the low energy theory, 
whereupon the SU(t) gauge theory becomes massive at low energies, with t vacua. From 
the M = 2 point of view, this mass generation in the SU(t) sector arises from t—1 condensed 
monopoles. To have t — 1 massless monopoles before the superpotential is turned on, the 
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right hand side of the equation of the TV = 2 curve y 2 = P t (x) 2 —4A 2t must admit a special 
factorization, 

P 2 (x) - 4A 2t = F 2 (x)H 2 _ 1 (x) 1 (2.42) 

where we can assume that Pt, F2, and Ht-i ah have leading coefficient 1. Moreover, setting 
F<i = (W') 2 + fo, we see that if W — x 2 /2, we want F 2 = x 2 — aA 2 for some constant a. This 
problem leads to t possibilities for P t and H t -i, found by using Chebyshev polynomials. 

Chebyshev polynomials of the first and second kind of degree t and t—1 respectively 
are defined as follows. (See Appendix C for some explicit formulas.) Set x = cos 9 and let 

/ x , „ . , , , . 1 d% , , sin(t6>) . d , 

%(x) = cos *0) U t -i(x) = --^(x) = —V' ( 2 - 43 ) 

t ax sine' 

the point being that cos(t6) and sin(t6*)/ sin 9 are both polynomials in x. Some identities 
that can easily be checked are 



T t (x) 2 -l = (x 2 -l)^_!(x) 2 
%(x) = \({x+ Vx 2 - If + (x - Vx 2 - 1)*) • 

From the second identityi it is clear that the coefficient of x l in T t (x) is 2 t ~ 1 . 
The solutions to ( ggg ) are! 

I ts-r- i X \ „ /x 9 .9.9 „ /x . , / x 



(2.44) 



P t {x) = 2K t rf% — - , F 2 (x) =x 2 -4^A 2 , H t _ x {x) =rf- L K t ~ l U, 



(2.45) 

with 77 a 2t-th root of unity, i.e. rj 2t = 1. Using the fact that Pt(—x) = (— l)*Pt(x), one 
can see that ( |2.45|) is a function only of r] 2 and hence gives us precisely the expected t 
solutions of the factorization problem. 

Using this, we can go back to the original problem of mapping U(N) solutions with 
a given superpotential, and some renormalization scale Aq, to U(tN) solutions with the 
same superpotential and renormalization scale A. 

By assumption, we have polynomials Pn(x), F2 n (x) and .£f/v_ n (x), all with leading 
coefficient 1, obeying the following relations: 

P 2 N (x) - 4A 2N = F 2n (x)H 2 N _ n (x) 

(2.46) 

F 2 n(x) = (W') 2 + f n - 1 (x). 



To prove this, write cos(t6) = \{e tei + e~ t9i ) 



2 

9 



We are making a slight refinement of the formulas as they have been presented previously. 
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Prom (|2.44j ), one can show that 2A tN r] t T t (x) with x = Pn( x ) satisfies 

2A tN rf%(x) = x tN + 0(x tN ~ 1 ) 
(2A tN V t T t (x)) 2 -4A 2tN = (P 2 (x)-A V 2 A 2N ) {r 1 t - 1 A N ^U t . 1 {x)^ . 

Using (|2.46| ) in the last equation with the identification Aq^ = r) 2 A 2N , one finds 



(2.47) 



(2A tN v t T t (x)) 2 - 4A 2 ^ = F 2n (x)H 2 N _ n (x) (y^ 1 A^^U^x)) 2 . (2.48) 

The factorization problem for vacua of the U(tN) problem with low energy gauge 
group U{l) n is 

P 2 N (x) - 4A 2tN = F 2n (x)H 2 N _ n (x). (2.49) 
Comparison of ( |2.48| ) and ( |2.49| ) leads to the identification of 

PtN(x) = 2A«Y* (^) 

F 2n (x) = F 2n (x) (2.50) 
H tN - n (x) = ^A^H^ixM^ 

as solutions of (|2.49| ). Since F 2n = F 2n , the vacua constructed this way for the U(tN) 
theory have the same superpotential as the vacua of the U(N) theory. 

For a given superpotential W, the U(N) theory has a finite number of vacua with 
given n. This number is independent of Aq. For t different values of the U(N) parameter 
Aq^, namely 

A 2N = V 2 A 2N , (2.51) 

we have been able to construct U(tN) vacua. Thus, this construction gives t times as 
many U(tN) vacua as there are U(N) vacua. The relation ( |2.51|) shows that, to obtain 



U(tN) vacua in which the 6*-angle equals 9tN, the analogous angle 9n of the U(N) theory 
must be djy = 6 t N + 2nl/t for some integer /. 

This construction gives t branches in U(tN) for every branch in U(N), since the con- 
struction depends on the choice of r] 2 . Since we will show that the construction multiplies 
the confinement index by t, the counting of branches is in accord with the observation in 
section 2.1 that Wilson and 't Hooft loops can distinguish t types of branch with confine- 
ment index t. 
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Unbroken Gauge Group 

As a check of our formalism, let us see how the tN vacua with unbroken U(tN) arise 
by applying the multiplication by t map to the unbroken vacua of U(N). 

In this case, Pn(x) satisfies (|2.46[ ) with n = 1, i.e. it is given by, 



P N (x) = 2A$e N T N 



x 



2eA 



with e 2N = 1. 

Using this in ( |2.50| ), we get 

Using Aq^ = n 2 A 2N , one can easily check that the combination squares to one and 

can be pulled out of % giving a factor of ( j . This implies that ( |2.52| ) can be written 

as 

P tN (x) = 2(eA ) tN T t (t n (^)) (2-53) 

Let us introduce a obeying a 2N = rj 2 as well as Ao = 3A, and define a = ea. Finally, 
using the identity %n(x) = %(Tn(x)) in ( [2.53| ), one gets 

P tN (x) = 2A tN a tN T tN (^L) . 

Recalling that r\ is a 2i-root of unity we get that a is a 2tiV-root of unity, and so is a since 
e is a 2iV-root of unity. This, as we wished to show, is the expected solution of the U{tN) 
problem. 

Generating Function of Observables Tr $ fc 

In section 2.1, a criterion for confinement was given in terms of the operator- valued 
one- form constructed from the generating function of Tr 

T(z) = ^Tr-^. (2.54) 

Let us denote by T t ]y(x) and T^{x) the one forms of U{tN) and U(N) respectively. 
Recall that (according to Appendix A of [fL7| ) 

Tr— L-\= - P ' n{x) (2.55) 
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From flQ3p and ( p30| ), we get 



P[ N (x) = 2 v t A tN ^T t '(x) = t^P'^A^-^Ut-iix). (2.56) 



Using this, ( gg7| ), and ( fT5q ), we get 



Tr -J— \ = t - P ^ (X) - (2.57) 
x - $ / 



'/*(*) -4A™ 

If we denote by $o the adjoint field in the U(N) theory, then ( |2.57|) can be written as 



x — <&/ \ x-$ 

or equivalently T t j^{x) = tT^{x). 

So the multiplication by t map multiplies all periods of T by t. Since the confinement 
index is the greatest common divisor of the periods, it follows that this map multiplies the 
confinement index by t. Since the Ni are certain periods of T, it follows also that the map 
multiplies the iVj by t, as originally shown in M. 

All The Confining Vacua 

As a special case, if we start with a Coulomb vacuum in U(N), the multiplication by 
t map produces vacua with confinement index t in U(tN). 

The converse is also true: all vacua in U(tN) with confinement index t arise in this 
way from Coulomb vacua in U(N). This can be shown by counting vacua. Because of the 
freedom to choose ry, the multiplication by t map produces t times as many vacua with 
confinement index t in U(tN) as there are Coulomb vacua in U(N). Let us show that this 
is the correct number. 

As explained in section 2.1, a Coulomb vacuum in U(N) with a classical limit based 
on Y\i U(Ni) is specified by picking some integers with < < Ni — 1 and such that 
the collection of the N and hi = — ri + i are relatively prime. Instead, to give a vacuum 
in U(tN) with classical limit based on Y\ i U(tNi) and with confinement index t, one must 
pick some integers r^, < fi < tNi — 1, such that the collection tNi and hi = r{ — r^i 
have greatest common divisor t. Any such ri are uniquely of the form = u + tr^ with 
< u < t — 1 (and that formula works for any u and for any that solve the problem in 
U(N)). As there are t choices of u, the number of U(tN) vacua with confinement index t 
should indeed be t times the corresponding number of U(N) Coulomb vacua. 
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2-4- Further Strong Coupling Generalities 

Here we will consider some additional issues from the strong coupling point of view, 
in which one starts with the solution of the M = 2 theory and treats the superpotential as 
a small perturbation. The SU(N) dynamics is described by the curve E: 

y 2 = P N (x) 2 - 4A 2N . (2.59) 

Here Pn(x) = det(x — $> c l)- Generically, the expectation value of $ breaks SU(N) to 
[/(l)^ -1 , and E has genus N — 1. Its Jacobian determines the effective couplings of the 
N — 1 low energy abelian vector multiplets. Including the U(l) center of U(N), which is 
free and not described by E, the low energy unbroken gauge group is 

If the right hand side of E has w double roots, this corresponds to existence of w 
independent massless monopoles. After small perturbation by a generic superpotential, 
reducing the supersymmetry to M = 1, the monopoles will condense, reducing the low 
energy gauge group from U(1) N to U(1) N ~ W . So to describe a vacuum of the Af = 1 
theory with low energy gauge group £7(l) n , we need N — n double zeroes of the polynomial 
P N (x) 2 -4A 2N . 

This polynomial factorizes 

P 2 N -4A 2N = (P N + 2A N )(P N -2A N ). (2.60) 

The two factors P N + 2A N and P N - 2A N never vanish simultaneously, so a double root 
is a double root of one factor or the other. To get w = N — n double roots, we therefore 
place double roots in one factor and s_ in the other factor, with w = s+ + s_. If w 
is given, there are w + 1 ways to pick s + and s_ with w = s + + s_, so if all possibilities 
actually occur, there are (at least) w + 1 branches of the moduli space of vacua on which 
the rank of the unbroken group is n = N — w. In the examples that we have looked at, 
all values of and s_ with 2s_|_,2s_ < iV are possible - they all are realized for some 
N th order polynomial Pn- (Note that every Pn(x) is of the form Pn(x) = det(x — $d) 
for some $ c ^.) We actually will meet cases with several branches for the same and s_. 

10 In this paper, we will not require $ to be traceless. So $ describes all chiral superfields, 
including the N = 2 partner of the center of U(N). But the Jacobian of S describes the gauge 
couplings only for the SU (N) gauge fields. When a superpotential is turned on, the trace part of 
<I> interacts with the other fields, but the U(l) gauge field remains free. 
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This multiplicity of branches cannot, as far as we know, be explained using conven- 
tional order parameters. The exchange s + <-> s_ can be understood as a discrete chiral 
symmetry of the Af = 2 theory under which changes sign. In some cases (with w 
close to N) that we will study in section 3, confinement goes a long way in explaining why 
there are several branches. However, in general, chiral symmetries and confinement (and 
any conventional order parameters we know) do not suffice to explain why there are so 
many branches. In section 4, we study examples in which w is small compared to iV and 
show that the branches with different s+ and S- can be characterized according to which 
holomorphic functions of the chiral order parameters vanish on a given branch. 

Now let us explain something that will be useful background for our computations in 
section 3. Consider a family of M = 2 vacua with N — n massless monopoles and given 
s+, S- (and s+ + S- = N — n). As we take A — > 0, can the theory converge to a classical 
limit in which the unbroken gauge group is Yl?=i U(Ni)? 

If such a limit does exist, then for very small A, the U(Ni) theories approximately 
decouple from each other. Each has its own polynomial Ppf. and its own mass scale Aj, and 
its own curve yf = Pjv 4 (#) 2 — 4A^ 7Vi . If each U(Ni) theory has Wi massless monopoles, we 
need Wi = N—n so that the full theory has N—n massless monopoles. But in U(Ni), the 
maximum possible number of massless monopoles is AT— 1, and since Xir=i(-^* — -0 = N—n, 
we must in fact have Wi = iVj — 1 for all i. 

The particular polynomials (x) such that Pjy. — 4A^ Ni has Ni — 1 double roots 



are known explicitly in terms of Chebyshev polynomials |2(J and were used in section 
2.3. Even without recalling the details, we can determine what numbers s+i and S-i of 
double roots of P/y. ± 2A? iVi are possible if s + i + s_j = N{ — 1. A polynomial of degree 
Ni has at most A^/2 double roots, so s+i, S-i < Ni/2. If Ni is odd, the only possibility is 
s+i = S-i = (Ni — l)/2, and if Ni is even, there are two cases (s+i, S-i) = (Ni/2, Ni/2 — 1) 
or (7Vy2-l,7Vy2). 

Finally we can answer our question of which classical limits can exist for given s + , s_. 
A branch of given s+ and s_ can have a classical limit with unbroken Y]a=i U(Ni) only if 
one can write s+ = J2i s - = J2i with < s+i, S-i < Ni/2 and s+i + S-i = ATj — 1. 
The above argument shows that this condition is necessary; in examples that we have 
studied, it is also sufficient. For example, in section 3, we will examine in detail the case 
of = 6 with n = 2, so s+ + S- = 4. For (s+,S-) = (3,1) or (1,3), our conditions 
only allow the classical limit U(4) x U(2), in agreement with what we will find. But for 
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(s + , s_) = (2, 2), they allow U(5) x U(l), U(4) x U(2), and C7(3) x U(3), all of which will 
duly appear. 

Intersections of Branches 

Now let us consider whether different branches might meet at singularities. For exam- 
ple, can a branch M. x with (s_|_, s_) = (a, 6) meet a branch A^2 with (s + , s_) = (o+l, 6—1)? 
This will naturally occur at a point with w + 1 massless monopoles and (s_|_, s_) = (a + 1, 6). 
Such a point is on both Aii and as well as being on a branch Ai with w + 1 massless 
monopoles. So in general, distinct branches with "adjacent" values of s+, s_ will meet at 
singular points with an extra massless monopole. 

Such a singularity corresponds to shrinking a cycle in the Riemann surface of the 
M = 2 theory. From ( |2.20| ) it is clear that the differential T cannot have a pole. Therefore, 
a cycle can shrink only if the period of T around it vanishes. This leads to restrictions 
on the possible singularities in a branch. Similar restrictions have been pointed out in || . 
They do not affect the examples we study below because even if § T around a certain cycle 
is not zero, we can find another cycle with vanishing period which can be shrunk. 

Here we have described the situation in an M = 2 language. From an M = 1 point of 
view, the situation is different. For a given branch with given w, the superpotential that 
leads to a given M = 2 curve on the branch can be found by the recipe of 0], reviewed 
in section 2.2. (It generically triggers condensation of all the massless monopoles.) For 
simplicity, we assume that on the branches M.\ and k = n so the superpotential is 
uniquely determined. The superpotential needed to get a given Af = 2 curve on a branch 
with (s+,s_) = (a, 6) is generically different from the superpotential needed to get the 
same curve on a branch with (s + ,s_) = (a + 1,6 — 1). So this type of intersection of 
branches is not really relevant for J\f = 1, which is our main interest in the present paper. 

What is more relevant for M = 1 is that a branch Aii with (s+, s_) = (a, 6) meets a 
branch M. with (s+, s_) = (a + 1, 6) at a point with an extra double root. (An example 
was discussed in 

The notion of M. and A^i "meeting" may seem confusing since from an M — 2 point 
of view, M. is just a subspace of M.\ singled out by existence of an extra double root. So 
let us elaborate. M. actually has a reduced rank of the low energy gauge group n = k — 1, 

11 The branch Jv[ likewise meets the branch .M2 with (s+,S-) = (a,b — 1); as we have just 
explained, this occurs with the same N = 2 curve but a different superpotential from the one that 
leads to Ali intersecting M. 
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so on this branch the superpotential is not uniquely determined by the Af — 2 curve, but 
depends on one extra complex parameter. This parameter controls the condensation of 
the massless monopole that appeared when Pjv + 2A 2Ar developed an extra double root. 
This parameter should be included in describing the A4 branch from an Af = 1 point of 
view. Thus, from an M = 1 point of view, A4 is not a subspace of .Mi, and in fact they 
have the same dimension. 

Starting on .Mi, the intersection with M. is obtained by adjusting the Af = 2 curve to 
get an extra double root. Starting from A4 (at a generic point where the extra monopole 
is condensed), the intersection is reached by varying the superpotential until it becomes 
the A4 i superpotential for the same Af = 2 curve. 

For this to make sense, we need to know that the A4\ superpotential of any Af = 2 
curve with (s + , s_) = (a + 1, b) is always one of the allowed superpotentials that can lead 
to that curve on AA. This can readily be shown from the generalized H recipe for the 
superpotential that was presented in section 2.2 (see ( |2.40|) and the discussion following 
it). ^ 

Sometimes it will happen that the branch A4 is confining while A4± is not confining 
(or has a smaller confinement index). If so, then as one approaches the intersection of the 
two branches starting on A4, the coefficient of the area law must vanish (as noted in an 
example in Q). 

Actually, a point with (s+, s_) = (a + l,b) is on an (o + l)-fold self-intersection of A4±, 
since there are a + 1 different ways to forget a double root of Pjy + 2A iV , leaving a of them. 
Such a point is likewise on a 6-fold self-intersection of AAi- However, the different branches 
do not really meet as Af = 1 theories, since they would require different superpotentials. 

Though in the present paper we focus on double roots of Pjy — 4A 2iV , it is also possible 
to consider the case of triple or higher order roots. The first example was studied by 
Argyres and Douglas |]21|| . At the Af = 2 level, before turning on a superpotential, these 
configurations are believed to correspond to non-trivial critical points. What happens if 
a superpotential is turned on, breaking to Af = 1? A generic superpotential will not be 
extremal at a point with a triple root, but by tuning parameters to get a superpotential 
that is extremal at such a point, one can presumably get a nontrivial critical point with 
Af = 1 supersymmetry. 



Off-shell Interpolation 
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Most of this paper is devoted to the question of branches of physical vacua. As the 
parameters in W($) are varied, the vacua change and we explore the various branches and 
the smooth interpolation from one vacuum to another within a branch. These are on-shell 
interpolations. It is interesting to ask whether there exists a Lagrangian which allows us 
to extrapolate between the different branches. 

This question can be answered at various levels. The most complete description of 
the theory is in terms of the microscopic gauge theory. It describes all the vacua of the 
theory. 

A more macroscopic description is the effective Lagrangian of [|J for the fields Si. The 
vacua are found by solving the equations of motion of Si. As we vary the parameters in 
W(<fr), these vacua change and we find the on-shell interpolation. One of our results is that 
we can continuously interpolate between vacua with one value of Ni to vacua with other 
values of Ni. Therefore, the superpotential of || has many stationary points corresponding 
to the different vacua with different N on the same branch. We will see examples of this 
in the next section. 

Another type of an effective Lagrangian we can consider is the following. The effective 
Lagrangian for a given branch has the form ( |2.13|) 



NiSi (log (J^j 2mbiSi + N *£>(SjSk) (2.61) 



(we absorbed tq into A and b\ = 0). We can replace it with 



J^Si (log ( ^£ ) +Ni ) +5^,0(5^*). (2.62) 



As we explained in section 2.1, this superpotential has more stationary points than ( |2.61| ) 



For small Si it has ]X Ni vacua, rather than a single vacuum. Therefore, it describes 
different branches, some of them can be confining and others can be in a Coulomb phase. 
A superpotential like ( |2.62f) might be a good effective Lagrangian for all vacua for a given 
N for small Si, i.e. at weak coupling. However, it is unlikely to describe the proper 
interpolation to large Si where vacua with different values of Ni are present. 
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3. Examples 

In this section, we will analyze several explicit examples demonstrating the phenomena 
we have discussed. In the first five subsections, we will study examples with gauge group 
U(N) for N = 2, . . . , 6 and with low rank of the low energy gauge group. The case that the 
rank is 1 is too trivial for our purposes; this means that U(N) is unbroken classically. The 
corresponding M = 2 curves are described by Chebyshev polynomials and were reviewed 
in section 2.3. We focus therefore on the next case that the low energy group has rank 
n = 2; thus, classically U(N) is broken to U(Ni) x U(N 2 ), and the question arises of 
whether smooth interpolation can occur between different pairs (Ni,N2). In section 3.6, 
we present some special examples with n > 2. 

For each N, the "new" vacua are the Coulomb vacua, since according to the discus- 
sion in section 2.3, the confining vacua with confinement index t are determined by what 
happened for U(N/t). In our systematic construction of all vacua with n = 2 for given 
N, we will come across all vacua, both Coulomb and confining. For the range of N we 
consider, interpolation between different pairs (Ni,N%) occurs only for Coulomb vacua, 
simply because (if t > 1) N/t is always too small for such interpolation to be possible in 
U(N/t). 

We will analyze the theory along the following steps. We start by considering the 
M = 2 U(N) gauge theory whose Coulomb branch is described by the hyper-elliptic curve 

y 2 = P N {x) 2 - AK 2N = (P N (x) + 2A N )(P N (x) - 2A N ) (3.1) 

where the N th order polynomial Pn( x ) = det(x — $ c z) parameterizes the point in the 
moduli space. This point is labeled by the N eigenvalues of the matrix $ c ; modulo per- 
mutations. 

To get n = 2, the N = 2 theory must have iV — 2 massless magnetic monopoles, so 
the polynomial ( |3.1| ) must have N — 2 double roots. Since P/v depends on iV complex 
parameters, the subspace on which P^ — 4A 2N has N — 2 double roots is two-dimensional. 
As discussed in section 2.4, the distribution of double roots between the two factors in (|3.1|) 
is labeled by two integers. There are s + double roots in the first factor and s_ double 
roots in the second factor. The two integers s+ and s_ must satisfy s + + s_ = N — 2 and 
s+, S- < N/2. Different values of s+ and s_ correspond to different branches. 
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After picking values of s + and s_ we solve the factorization problem 

P N (x) + 2A N = H 2 (x)R N - 2s+ {x) 

~ ~ (3-2) 
P N (x) - 2A N = H 2 (x)R N - 2s _ (x) 

with polynomials H s+ , H s _, Rn- 2s+ and Rn- 2s _, where the subscript of the polynomial 
denotes its degree. We normalize the coefficient of the highest power of x in each polynomial 
to be one, and solve for the other coefficients such that 

H 2 s+ (x)R N - 2s+ (x)-4A N = H 2 _(x)R N - 2s „(x). (3.3) 

We will see in the examples below that often there are several distinct branches for the 
same (s_|_, s_). 

The moduli space M. is parametrized by the coefficients in the polynomials H s+ , 
H s ^, Rn-2s + and Rn- 2s . As we mentioned above, M. is two dimensional. In solving the 
factorization problem ( |3.2j ), we will sometimes shift a; by a constant in order to simplify the 
equations. This constant can trivially be reinstated at the end of the calculation. Ignoring 
this constant, our moduli space is one dimensional. 

We then examine the semiclassical limit A — ► 0. Since we limit ourselves to the case 
n = 2, in the classical limit the microscopic U(N) gauge group is broken to U(Ni) x U(N 2 ). 
We want to determine for each branch M. a of the moduli space what values of N\ and N 2 
are possible. 

We will also study singular points on M. a . Triple roots of y 2 lead to Argyres-Douglas 
points [^TJ. We will not determine them here. We will focus on singularities in which y 2 
has a new double root. These occur when either Rn_ 2s (x) or Rj^_ 2t {x) acquire a double 
root. This means that these points have (s_|_ + 1, s_) or (s + , s_ + 1) double roots in the 
two branches. At these subspaces, there are — 1 massless monopoles. As explained in 
section 2.4, at these points M. a meets a branch with n = 1, i.e., the fully confined branch 
that classically has unbroken U(N). 

We then break M = 2 to M = 1 by turning on a tree level superpotential which leads 
to a vacuum at a point on A4 a ; this lifts the degeneracy of -M a , and only a finite set of 
vacua survive for given superpotential. We want to determine what superpotential can 
lead to a given point on the moduli space. Since n = 2, the superpotential must have 
at least two critical points, and to keep things simple we will consider the case of a cubic 
superpotential, which has exactly two critical points. In section 2.2, we reviewed the recipe 
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of H for determining what cubic superpotential, normalized to have leading term x 3 /3, 
leads to a given point on the M = 2 moduli space. We take the matrix model curve, which 
in this case is 

Vm = Rn-2s+(x)Rn-2s_(x), (3.4) 



and we write 



R N -2s + (x)R N -2s. (x) = W'{xf + h{x). (3.5) 



Recall that + = iV — 2, and therefore the left hand side of ( |3.5|) is a quartic polyno- 
mial. ( |3.5|) determines W', and hence determines W up to an irrelevant additive constant. 
Moreover, the value of the gluino bilinear S can be read off from fi(x); if fi(x) = Ax + B, 
then S = Si + S 2 = -A/4. 

The intersection points of different branches are points where has a double root 
and therefore equation ( |3.4|) describes a genus zero curve. This is consistent with the fact 
that such points are also on the n = 1 branch because, as explained in section 2.2, when 
n < k (here n = 1 < k = 2), the matrix model curve has double zeros. The physics of 
these points in the Af = 1 theory is interesting. The low energy spectrum at these points 
consists of the two U(l) multiplets which exist at generic points in M. a as well as another 
massless monopole. Moving away from this point along M. a , this monopole acquires a 
mass. Moving away from this point through the condensation of this monopole takes us to 
the fully confining vacua with unbroken U(N) and only a single massless U(l) multiplet. 

If we denote the critical points of W as x\ and x%, we have 

W' = {x-x 1 ){x-x 2 ). (3.6) 

The physics is essentially unchanged by adding a constant to x, which shifts x\ and x 2 by 
a constant. The invariant information contained in the choice of W , modulo shifting x 
and exchanging the two critical points, is contained in the parameter 

A = (xi-x 2 ) 2 . (3.7) 

It is possible and sometimes convenient to shift x so that a superpotential of given A takes 
the form 

x 3 A 

W=---x. (3.8) 

W is odd under the transformation x —>■ —x, which thus corresponds to an i?-symmetry. 
(Of course, the symmetry is present whether or not we shift x to put W in the form (|3.8| ).) 
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The critical points xi, X2 are exchanged by the symmetry, so a choice of classical vacuum 
with Ni 7^ N 2 spontaneously breaks the symmetry. In some examples, the symmetry 
exchanges different branches Ai a , and in other cases, it acts within a fixed branch. 

The matrix model formulas can be used to compute T(x). Though we will not do this 



in detail, it could be done as follows. Following Jr? 



T (x) = /tt -^-\ = - P ' n(x) -.dx. (3.9) 

From O 

P' N (x) =H s+ (x) (2H> s+ (x)R N - 2s+ (x)+H s+ (x)R' N _ 2s+ (x)) 

= H s _ (x) (2H' S _ (x)R N _ 2s _ (x) + H s _ {x)R' N _ 2s _ (x)) (3.10) 
= -^c(x)H s+ (x)H s _(x). 

Since we learn from the first two expressions that P' N (x) is divisible by H s+ (x) and H s _ (x), 
c(x) must be a polynomial (the factor of —1/4 is in order to agree with the conventions of 
|fP7|1 ). Therefore ( |3.9| ) can be written as 

T(x) = -^-dx (3.11) 

which leads to the identification c(x) = -4 [W'(x)T(x)] + = -4Tr i n p. 

In some of our examples, we will find S = 0, and find that this can be interpreted as a 
consequence of the symmetry. However, the symmetry is special to a cubic superpotential. 
A generic superpotential of higher degree does not have such a symmetry and it can still 
lead on certain branches to S = 0, as we will demonstrate in section 3.6. Lacking such a 
symmetry, vanishing of S on a branch of vacua would not be natural in a nonsupersym- 
metric field theory. With supersymmetry, such "unnatural" phenomena can be dictated 
by holomorphy. In section 4, we discuss more elaborate examples of such phenomena. 

Let us consider the examples. 

3.1. U{2) 

Our first example is a U (2) gauge theory. The hyper-elliptic curve (|3.1|) is 



y 2 = (P 2 (x) + 2A 2 )(P 2 (x) - 2A 2 ). (3.12) 
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The factorization problem (3^2) in this case is quite trivial, since to get n = 2 we do not 
need double roots at all - just the classical symmetry breaking U(2) — > U(l) X U(l). We 
write 

P 2 (x) = (x-a)(x-b) (3.13) 

Our moduli space A4 with n = 2 is parametrized by a and b modulo the exchange a <-> 6. 
In other words, it is the whole M = 2 Coulomb branch. 

The semiclassical limit A — > is also trivial since P2(%) = (x — a)(x — b) is independent 

of A. Setting P 2 (x) =det(x-$ c i), we find ($ cI ) = ^ °Y breaking 17(2) to C/(l)xC/(l). 

The subspace with a massless monopole is determined by looking for points where y 2 
has a double roots. This happens when 

(a-6) 2 = ±8A 2 . (3.14) 

These subspaces are also on the unbroken £7(2) branch corresponding to n = 1. If we shift 
x such that a = —b (as in SU(2)), at these points 

P 2 (x) = x 2 ^2A 2 , =a; 2 (x 2 T4A 4 ). (3.15) 

Note that /^(a?) = x 2 — 2p 2 A 2 with p 4 = 1; this is equal to 2K 2 p 2 T 2 {^^) from section 
2.3. After a small perturbation causing the monopoles to condense, one of these points 
has confinement and one has oblique confinement. 

We now break Af = 2 to Af = 1 by turning on a cubic superpotential chosen to put 
the system at the point (a, b). We find it through the matrix model curve 

yl = P 2 - 4A 4 = ((x -a)(x- b)f - 4A 4 (3.16) 

from which we derive 

W'(x) = (x-a)(x-b) 
f{x) = -4A 4 (3.17) 
S = 0. 

In this case, vanishing of S can be attributed to the discrete chiral symmetry. 

If W is given, ( |3.17| ) determines a and b modulo the exchange a <-> b. Since the 
moduli space is obtained by dividing out by this exchange, there is in this example only 
one vacuum for given W. This occurs whenever all Ni are 1, as then there is no strong 
dynamics at low energies. 
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3.2. U{3) 

Our next example is somewhat more interesting. It is based on the gauge group £7(3). 
We start by finding the subspace of the M = 2 moduli space with one monopole. The 
Z7" (3) curve 

y 2 = (P 3 (x) + 2A 3 )(P 3 (x) - 2A 3 ) (3.18) 

should have a single double root. The solution of this problem has two branches labeled 
by rj = ±1 depending on which factor of ( |3.18j ) has the double root; i.e. (s+, s_) = (0, 1) 
or (1,0). In either case, the factorization we need is 

P 3 (x) + 2r)A 3 = (x-a) 2 (x-b), (3.19) 

implying also 

P 3 (x) - 2r/A 3 = (x- af{x - 6) - 4r/A 3 . (3.20) 

As in the previous example, a and b are the parameters on the moduli space; now there is 
no symmetry exchanging them. 

The semiclassical limit A — > leads to P 3 (x) — > (x — a) 2 (x — b) and therefore it 
corresponds to £7(3) broken to U{2) x £7(1). The existence of two branches - via the choice 
of 77 - has a simple explanation: the SU(2) that is unbroken classically is strongly coupled 
at low energies and has two supersymmetric vacua, corresponding in our problem to the 
two branches. 

Each branch has three one-dimensional subspaces with one more massless monopole. 
They are determined by looking for points where y 2 has two double roots. This happens 
when 

(a-b) 3 = 27r/A 3 . (3.21) 

Since these points have (s+,s_) = (1,1), at these points the two branches labeled by r\ 
meet each other, and they also meet the n = 1 subspace of the Af = 2 moduli space. 
Therefore, they can be interpreted as unbroken £7(3) vacua, in keeping with the discussion 
in section 2.4. In fact, once we impose ( |3.21| ) in the form (a — b) = 3pA with p 6 = 1 and 
shift x — ► x + I (2a + 6), P3 becomes a; 3 — 3p 2 A 2 x which is equal to 2A 3 p 3 73( 7^) from 
section 2.3. 

We now break M = 2 to M = 1 by turning on a superpotential which puts the system 
at the point (a, b) in the Af = 2 moduli space. We find it through the matrix model curve, 
as in (13.51): 



y 2 m = {x- b) ({x - a) 2 (x - b) - 4r/A 3 ) = ((x - a)(x - b)f - 4r/A 3 (x - b). (3.22) 
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From this we derive 

W'(x) = (x-a)(x-b) 
f(x) = -4nA 3 (x-b) (3.23) 
S = V A\ 

For fixed superpotential, there are two values of (a, b). Either of them can be the 
double root of P3 + 2??A 2 , and the vacuum also depends on the choice of n. So for each tree 
level superpotential, the system has four vacua with low energy gauge group U(l) 2 . Semi- 
classically, they can be interpreted as follows. First, the unbroken £7(2) can be associated 
with either one of the two different minima of the potential, spontaneously breaking the 
global Z 2 symmetry. Each of these choices leads from the strong coupling SU(2) dynamics 
to two vacua, which differ by the sign of S = nA 3 . 

3.3. U{4) 

We now repeat the same analysis for U(4). We look for the subspace of the Af = 2 
theory with TV — n = 4 — 2 = 2 magnetic monopoles. The two double roots of the M = 2 
curve 

y 2 = (p 4 ( x f + 2A 4 )(P 4 (x) 2 - 2A 4 ) (3.24) 
can be distributed between the two factors as (s + , s_) = (1, 1), (2, 0) or (0, 2). 

Confining Branches, Monopoles Distributed As (s + , s_) = (2,0) or (0,2) 

With both double roots in the same factor of Pf — 4A 8 , we get two branches, labeled 
by the choice of n = ±1: 

P 4 (x) + 2nA 4 = (x 2 -ax + b) 2 

(3.25) 

P 4 (x) - 2nA 4 = (x 2 - ax + b) 2 - 4??A 4 . 

The semiclassical limit of A — > leads to P 4 (x) — > (x 2 — ax + 6) 2 , which has two 
double roots, so U(4) is broken to U{2) x U{2). That this is the only classical limit for 
(s + , S-) = (2, 0) or (0, 2) could be predicted from the reasoning in section 2.4. 

Each branch has four singular subspaces with another massless monopole, i.e., on 
which P 4 (x) — 2rjA 4: also has a double root. This happens at the solutions of 



(46 -a 2 ) 2 = 64 V A 4 . 
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(3.26) 



These points correspond to (s + , s_) = (2, 1) or (1, 2), so they are on the branch with n = 1 
and unbroken £7(4). 

As in (13.51), the matrix model curve is 



Vr. 



(x 2 - ax + b) 2 - 4rjA 4 



(3.27) 



from which we derive 



W'(x) = x 2 — ax + b 
f(x) = -At]A 4 
s = 0. 



(3.28) 



We see that for each superpotential, we have two vacua (labeled by rj), one on each branch. 
This can be understood semiclassically as follows. The low energy U(2) x U(2) is charac- 
terized by scales which are given in the semiclassical limit by Af = A 2 ~ a ^-4f Gluino 
condensation in these groups leads to four vacua Si ~ ±Af , S2 ~ ±Af . The Z2 exchanges 
Si and S2, and since it is an .R-symmetry, it acts by Si <-> — S2. There should therefore be 
two Z2-invariant vacua, with Si = — S2 and S = 0; these are evidently the vacua that we 
have found on the branches with (s+, s_) = (2, 0) or (0, 2). We will clearly have to find 
the two vacua with Si = +S2 and S 7^ elsewhere. 

Of the four U(2) x U(2) vacua, the analysis of section 2.1 shows that two are confining 
and two are in a Coulomb phase. In fact, the vacua that we have found are confining, 
because the branch that we have just described can be constructed by "multiplication by 
2," of P2(x) = x 2 — ax + b as in section 2.3, 



Pt(x) = 2 A V T 2 



(3.29) 



with p 4 = 1 and p 2 = rj. This is the same as ( 3-25 ). 



Coulomb Branch, Monopoles Distributed as (s+,s_) = (1,1) 

The other branch occurs when each of the factors P±{x) ± 2A 4 has a single double 
root. It is not too hard to show that, modulo the freedom to add a constant to x, this 
implies that 



[x — a) 



P 4 (x) + 2 A 4 
P 4 (x) -2 A 4 = (x + a) 



4 



[x + a) 2 + ^ T (x + 2a) 
a" 3 

(x — a) 2 H t(x — 2a) 

a 6 



(3.30) 
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0: 



[x - a) 2 (x + a) 2 . Therefore £7(4) is broken to £7(2) x U(2). 



for some a. 

There now are two semiclassical limits with A 

1. Fixed a: Pa{x) 

2. A, a — > with fixed v = A 4 /a 3 : P^x) — > x s (x + v), showing that here £7(4) is broken 
to £7(3) x U(l). 

This is the first example of our new duality: the same moduli space has two distinct 
semiclassical limits corresponding to different gauge groups. One can continuously deform 
a £7(2) x £7(2) vacuum to a £7(3) x £7(1) vacuum. 

Since the numbers 1 and 3 are relatively prime, these £7(3) x £7(1) vacua are all in 
a Coulomb phase. This whole branch is therefore in a Coulomb phase, and the classical 
limit with £7(2) x £7(2) will turn out to give the two £7(2) x £7(2) vacua that did not appear 
earlier. The upshot will be that £7(2) x £7(2) vacua that are in a Coulomb phase can be 
smoothly transformed to £7(3) x £7(1) vacua, while confining £7(2) x £7(2) vacua cannot 
make such a transformation, since there are no confining £7(3) x £7(1) vacua. 

The special points on this branch where another monopole becomes massless are at 
the solutions of 

16a 8 = A 8 . (3.31) 

Explicitly, a = ^pA with p 8 = 1. Note that shifting x — > x — and using the solutions 
to ( p.31|) we get from (|3.30| ) that P±(x) = 2A 4 p 4 74( t^-). These are the same four points 
that we found in ( |3.26|) . 

From ( |3.3U| ), we find the matrix model curve, as in ( |3.5| ): 



V 2 m 



A 4 

(x + a) 2 H t{x + 2a) 



a. 



(x 



x 2 + x 



A 4 



a 



4A 4 



,2 A 4 , 
-a) H — 5- (x — 2a) 

a 6 

4A 8 



(3.32) 



-x 



The superpotential and f(x) are therefore 



A 4 



W'ix) = x + x—^ - a 



4A 4 4A 8 

X r- 



a, 



(3.33) 



a 
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Finally, we count the number of vacua for fixed tree level superpotential. The most 
convenient way to proceed is to compute the discriminant of the quadratic polynomial 
W' = x 2 + x(A 4 /a 3 ) - a 2 ; it is 

A 8 

A = 4a 2 + — . (3.34) 
a 

By shifting x, we could put the superpotential in the form (|3~8|). If A is given, then ( |3.34| ) 
is equivalent to an eighth order equation for a, namely 

4a 8 - Aa 6 + A 8 = 0. (3.35) 

So there are eight possible values of a for every A, implying that the theory has eight 
vacua for each tree level superpotential. In the semiclassical limit of large A/ A 2 , these 
vacua can be understood as follows: 

1. U(3) x U(l) leads to 3 x 2 vacua. The factor of 2 arises from the broken global 
symmetry which exchanges the two minima of the potential. Indeed, for large A/A 2 , 
( EP51 ) has six roots with a ~ (A 8 /A) 1 / 6 . 

2. U(2) x U(2) leads to 2x2 — 2 = 2 vacua. Here we do not need to include the factor of 
2 which exchanges the two minima of the potential and we subtract the two confining 
vacua that we found on other branches. To check this prediction, we see that for large 



A/ A 2 , ( EP5D has two roots with a ~ (A/4) 



1/2 



3.4. U(5) 

To find points on the Coulomb branch for U(5) with three massless monopoles, we 
want polynomials P5 such that 

y 2 = (P 5 (x) + 2A 5 )(P 5 (x) - 2A 5 ) (3.36) 

has A^ — n = 5 — 2 = 3 double zeros. This leads to two branches with (s+, s_) = (2, 1) or 
(1, 2). We will label them by rj = ±1. P5 + 2r]A 5 should have a pair of double roots, and 
P5 — 2?yA 5 should have a single double root; we use the freedom of shifting x to place this 
double root at the origin. Requiring the existence of the stated double roots leads to 

P 5 (x) + 277A 5 = (x 2 + ax - 2ac) 2 (x + c) 

(3.37) 

P 5 (x) - 2r]A 5 = x 2 [x 3 + (2a + c)x 2 + a(a - 2c)x - ac(3a + Ac) 



where the parameters a and c are constrained to satisfy 

a 2 c 3 = rjA 5 . (3.38) 
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From ( 3.38 ) it is clear that there are two classical limits: as A — > 0, either o -> or 
c — > 0. In the former case, P§(x) — > x 4 (x + c) corresponding to U(4) x £7(1) while in the 
latter, Ps(x) — > x 3 (x + a) 2 giving U(3) x U(2). This is another example showing that the 
same branch can have different classical limits. (If we let both a and c to to zero, we will 
get Ps(x) = x 5 , corresponding to a classical limit with unbroken Z7 (5).) 

The points with another massless monopole correspond to (s+,s_) = (2,2). They 
occur when P*>(x) — 2r]A 5 has two double roots. This happens when a 2 + llac — c 2 = 
subject to the constraint ( |3.38| ). Solving for c we get, 



c=-(ll + 5eVb)a with e 2 = 1 (3.39) 
2 



and 

a 



5 - (-682 + 305V5e)?7A 5 = (-2 + eV^A 5 . (3.40) 



Equation (|3.40| ) has 10 solutions for each 77 

a = ^(-2 + ev / 5)A, £ 5 = 77. (3.41) 

However, using the freedom to shift x such as to cancel the x 4 term in Ps(x) (bring it to 
the SU(5) form), x = x — |a(3 + \/5e), (|3.37| ) becomes 



P 5 (x) + 2ryA 5 = (x + 2£A) (x 2 - ^Ax - £ 2 A 2 ) 

o (3-42) 

P 5 (x) - 277A 5 = (x - 2£A) (x 2 + £Ax - e 2 A 2 ) 

from which it is clear that there are only 5 different points on each branch. Note that this 
is equal to P${x) = 2A 5 ^ 5 T 5 (^) from section 2.3. 

Why has the same P5 appeared for two different values of the pair (a, c)? These values 
of (a, c) are characterized by the fact that Ps(x) — 2r)A 5 has a second double root - but 
there are two choices of which double root is the "second" one. 

Using the recipe of (|3.5[) , we identify the matrix model curve, 

y 2 a = (x + c) [x 3 + (2a + c)x 2 + a(a - 2c)x - ac(3a + 4c)] . (3.43) 
From this we find 

W'(x) = x 2 + (a + c)x — ac 
f(x) = -4c 2 a(x + a + c) (3.44) 
S —c 2 a. 
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Finally, we count the number of vacua for fixed tree level superpotential, by the 
same reasoning as in previous examples. By computing the discriminant of W , we get 
A = a 2 + 6ac + c 2 . Using ( |3.38| ) to solve for c, we find that for fixed A there are 20 vacua, 
corresponding to the roots of 

a 8 (a 2 - A) 6 - 2a 4 (17a 2 + A) (1153a 4 + 142a 2 A + A 2 ) A 10 + A 20 = 0. (3.45) 

In the semi classical limit of large A/ A 2 , these vacua can be understood as follows: 

1. 17(4) x £7(1) leads to 4 x 2 = 8 vacua, with a ~ exp(±27rz/8)(A 10 /A 3 ) 1 / 4 . 

2. £7(3) x £7(2) leads to 3 x 2 x 2 = 12 vacua, which are at roughly a ~ ±\/A. 

We have not found any confining branches for £7(5) for the simple reason that 5 is a 
prime number, so if 5 = N\ + A^, then N\ and A^ are relatively prime. Hence for £7(5), 
all vacua are in a Coulomb phase. 

3.5. £7(6) 

£7(6) is the last and richest example that we will examine in detail. We look for points 
in the Coulomb branch of the M = 2 theory with N — n = 6 — 2 = 4 double zeros. These 
can be distributed between the two factors of 

y 2 = P 6 (x) 2 - 4A 12 = (P 6 (x) - 2A 6 )(P 6 (x) + 2A 6 ) (3.46) 

as = (3,1), (1,3) or (2,2). 

Confining £7(2) x £7(4) Branches, Monopoles Distributed As (s_|_, s_) = (1,3) or (3, 1) 

We first wish to consider the branches of the moduli space of £7(6) with four massless 
magnetic monopoles distributed as (3, 1) and (1,3) between the two factors. The solution 
of the factorization problem (|3.2|) is 



P 6 (x) + 2nA 6 = \(x - a) 2 (x - b) - 2eA 3 l 

(3.47) 

P 6 (x) - 277A 6 = (x- a) 2 (x - b) [(x - a) 2 (x - b) - 4eA 3 ] 

with n 2 = 1 and e 2 = n. 

The semiclassical limit A — > leads to Pq{x) — > (x — a) 4 (x — 6) 2 , and therefore £7(6) 
is broken to £7(4) x £7(2). That this is the only semi-classical limit with (s+, s_) = (3, 1) 
or (1, 3) follows from the reasoning in section 2.4. 
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The points with another massless monopole occur when Pe(x) — 2r/A 6 has two double 
roots. This happens when (a — b) 6 = (3A) 6 . These special points on the (s+, s_) = (3, 1) 
branches have (s+,s_) = (3,2); they also meet the (s+, s_) = (2,2) branches we will 
discuss below and the U(6) branch with n = 1. The special points on the (s+, s_) = (1, 3) 
branches have (s+, s_) = (2, 3) and they meet the (s+, s_) = (2, 2) branches and the n = 1 
branch. 

Using ( |3.5[ ), the matrix model curve is found from ( |3.47| ) to be 



(x-af(x-b) 2 -4eA 3 (x-b). (3.48) 



From this we derive 



W'{x) = (x-a)(x-b) 
f(x) = -4eA 3 (x-6) (3.49) 
S = eA 3 . 

For fixed superpotential A = (a — b) 2 there are 8 solutions. They correspond to the 
four solutions of e 4 = 1 for each choice of (a — b) = ±yA. In the semiclassical limit 
these 8 vacua can be understood as follows. U(4) x U{2) gives a total of (4 x 2) x 2 = 16 
vacua, where 4x2 reflects the strong dynamics and a factor of 2 comes from the broken 
7i2 symmetry. But, following the analysis in section 2.1, only half of these vacua (those 
with b even) are confining. So we expect 8 confining vacua. 

The eight vacua with (s+,s_) = (3,1) or (1,3) actually are confining, since they 
can be obtained by applying the "multiplication by 2" map to the polynomial Ps(x) = 
(x — a) 2 (x — b) — 2pAj] with p 2 = 1, which describes the breaking U(3) — > U (2) x U (1) (and 
was presented in ( |3. 19|) ) . Following the general discussion of section 2.3, the scales of the 
theories are related by A(j = cr 2 A 6 with a 4 = 1. Using this in P 6 (x) = 2A 6 a 2 T 2 (^^) we 
recover ( [3.471) with e = pa. 



Confining U(3) x U(3) Branches, Monopoles Distributed as (s_|_,s_) = (2,2) 

We now consider the case of (s+, s_) = (2, 2). In this case, the factorization problem 
( |3.2| ) has two types of solution. The first case leads to three branches parametrized by a 
cube root of unity p 

P 6 (x) + 2A 6 = (x 2 + g - pA 2 ) 2 (x 2 + g + 2pA 2 ) 

P 6 (x) - 2A 6 = (x 2 + g + pA 2 ) 2 (x 2 +g- 2pA 2 ) (3.50) 

P 6 (x) = (x 2 + s)[(£ 2 + s) 2 -3p 2 A 4 ]. 
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These three branches are distinct. Since Pq depends on A 4 and not on A 2 , the change 
A 2 — > —A 2 (and therefore A 6 — > —A 6 ) does not lead to more branches. Instead, there is a 
global symmetry g — > —g which acts on each branch. 

In the classical limit A — > we have Pe(x) — > (x 2 + g) 3 and therefore {7(6) is broken 
to {7(3) x {7(3). On each of these branches there are two "£7(6) points" - points with five 
double roots: g = ±2pA 2 . 

We now break M = 2 to M = 1 by turning on a superpotential W which puts the 
system in the ground state labeled by g. It is determined from the matrix model curve, 
which is found, as in (|3.5| ), from (ft.50|) : 



y 2 m = (x 2 + g + 2pA 2 )(x 2 + g- 2pA 2 ) = (x 2 + gf - 4p 2 A 4 . (3.51) 
From this, we find 

W'ix) = x 2 + g 
f(x)= -4p 2 A 4 (3.52) 
S = 0. 

These expressions are consistent with the Z2 R-symmetry which maps $ — > — $ and 
S -> -S. 

From the expression W = x 2 + g, it is clear that W determines g and hence Pn, 
so for fixed W, there is a single quantum vacuum on each branch. The three branches 
give a total of three vacua. Semiclassically, these can be interpreted as coming from 
L^(6) — > U (3) x i7(3), with the two gluino condensates in the two SU (3) factors anti-aligned. 
As in the discussion after ( [Tgp of {7(4) -> {7(2) x {7(2), this leads to S = Si + S 2 = 0, 



which respects the global Z2 R-symmetry of the system. 

{7(3) x {7(3) is expected to give a total of 3 x 3 = 9 vacua, but according to the analysis 
in section 2.1, only 3 of them (those with 6 = 0) are confining. We claim that the three 
vacua obtained by the construction above are confining. In fact, they can be obtained by 
the "multiplication by 3" of i^^) = x 2 + g which is a solution to {7(2) — > {7(1) x {7(1). 
One can check that P 6 (x) = 2A 6 e 3 T 3 (§^) reproduces flOcD with e 6 = 1 and e 2 = p 2 
cubic roots of unity. 

The fact that these three branches are confining, with a confinement index of 3, 
explains why they have no semi-classical limit with {7(4) x {7(2) or {7(5) x {7(1). There 
are three branches for a reason explained in section 2.1; they are distinguished by which 
of W U H, for u = 0, 1, 2, has no area law. 
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Coulomb Branch, Monopoles Distributed as (s+,s_) = (2,2) 

The second kind of solution of the factorization problem Q3.2|) in this case is 



P 6 (x) +2A e 



x 2 + (ft + g)x + 



(3ft + g)(9h 3 + 15h 2 g - hg 2 + g 3 )^ 2 
108ft 2 



x 



(h-g)(3h-g) 2 {Sh + g) 



108ft 2 



P fi (x)-2A C 



2g 



& + -£-r + (h- g)(x+ -§-) + 



2g 2 (h + g)(3h + g) 2 (3h-g) 
[X+ 3 } 108ft 2 



2g {3h- g){9h 3 -lhh 2 g -hg 2 - g 3 



108ft 2 



with gf and ft satisfying the constraint 



(3.53) 



g 5 (g 2 - 9h 2 ) 2 = 27 3 ft 3 A 6 . 



(3.54) 



We used a parametrization which makes the global symmetry g — > —9 combined with 
A 6 — ► —A 6 manifest. 

The classical limit A — > is obtained at: 

1. g -»• with finite ft. Here P 6 (x) -> (x + §) 5 (x - J); i.e. £7(6) -> £7(5) x £7(1). 

2. # -> 3ft with finite ft. Here P 6 (x) -> x 2 (x + 2ft) 4 ; i.e. £7(6) -> £7(4) x 17(2). 

3. g -> -3ft with finite ft. Here P 6 (x) -> (x - 2ft) 2 x 4 ; i.e. 17(6) -> Z7(4) x £7(2). 

4. ft, </ -> with v = g 2 /h fixed. Here P 6 (x) -> (x 2 + f^) 3 ; i.e. (7(6) -> U(3) x (7(3). 
This is a richer example demonstrating that the same branch can have different clas- 
sical limits. Since 1 and 5 are relatively prime, what we have found is clearly a Coulomb 
branch. 

The £7(6) points are obtained when either Pq(x) +2A 6 or Pq(x) — 2A 6 acquires another 
double zero. This happens when 



(h + g)(3h + g) 2 (3h-g) 
108ft 2 







(3.55) 



or the same equation with g — > —g. Taking into account (|3.54j ), this happens only for 



3a6 



ft = g with g a solution of (2g) b = 27 3 A 



(3.56) 



or 



ft = —g with g a solution of (2gY 



3a6 



-27^A 



(3.57) 
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Superficially, this leads to 12 points. However (after an appropriate shift of x) one finds 
that at these points Pq(x) ± 2A 6 depend only on g 2 and therefore there are only 6 such 
points. Indeed, both solutions lead to P§{x) = 2 A 6 p 6 %(-^x) w ^ n P 12 = 1- At these 
points, this branch meets the various confining branches that we have found. 

We now break Af = 2 to Af = 1 by turning on a superpotential W. It is determined 
from the matrix model curve 



Via 



■ (h-g)(3h-g) 2 (3h + g) 
108/i 2 



( | 2g 2 (h + g)(3h + g) 2 (3h-g) 
[X + 3 } 108/i 2 



2 (7 4 -6(7 2 /i 2 -27/i 4 V 4g 2 (g 2 -9h 2 ) 2g 2 (g - 3h) 3 (g + 3h) 



x2+ 3 9x+- J W8h2 J " " glh 'x+ 729h2 

(3.58) 

and therefore 

2 2g g* - 6g 2 h 2 - 27/i 4 „ 2 fir 4 - 18^ 2 /i 2 - 27/i 4 

(X) = X + —X H — — — tt; = X + 



3 108/i 2 108fr 2 

% 2 (^ 2 -9/i 2 ) , 2^ 2 (^-3/i) 3 (^ + 3/i) V(^ 2 -9/i 2 )^, 2^ 2 (^ 4 -81/i 4 ) 



= oTT x + ^3. = oTT x + 

S 



81h 729h 2 81h 729/i 2 

g\g 2 - 9h 2 ) 



81h 

(3.59) 

We used the freedom to shift x by x = x — % to put the superpotential W in a canonical 



3 

form W = ^- — ^x. 

We would like to parametrize the theory in terms of A. Including the constraint 
( |3.54|) , we have two equations with two unknowns: 

_ 27/i 4 - g 4 + 18g 2 h 2 

~ 27ft 2 (3.60) 

g 5 (g 2 - 9h 2 ) 2 = 27 3 /i 3 A 6 . 



Using these equations we can solve for h, 



V(i6^Q(^ 2 -|A) 2 -3(|A)(3A) 12 ) 
9(3A) 6 (32(7 10 ((7 2 - |A) +3(3A) 12 ) 1 ' ' 



and find a single equation for a single unknown g, 



1» ( 9 2 - ?aV = (3A)« (Is 12 - | (§a)V + i ( V) V + JL(3A)» | . (3.62) 



In deriving this equation we assumed that g ^ 0. 
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Since the final constraint ( j3.62|) is of degree 24, 



for fixed superpotential (fixed A) there 



are 24 vacua. These can be interpreted as follows: 

1. U(l) x U(5) vacua. Their number is 1 x 5 x 2 = 10 (the factor of 2 is associated with 
the global Z 2 symmetry of the theory). 

2. U(2) x U(4) vacua. Their number is 2x4x2 — 8 = 8 (we subtracted the 8 confining 
vacua) . 

3. U(3) x U(3) vacua. Their number is 3x3 — 3 = 6 (we subtracted the 3 confining 
vacua, and we did not multiply by 2 because for U(3) x U(3) the Z 2 symmetry is not 
broken) . 

It is interesting to contrast this situation with the classical theory (A — > 0). Here we 
find 3 solutions: g = 0, ±|\/A- The corresponding values of h depend on the details of 
the limit. The two solutions with g = ilv^A correspond to U{2) x £7(4) vacua (compare 
with the analysis of the semiclassical limits above). Each of them splits to 4 vacua in the 
quantum theory. The solution (7 = corresponds to U(l) x U(5) and to U(3) x U(3). 
This point splits in the quantum theory to 16 vacua (compare with the analysis of the 
semiclassical limits above). 

3.6. Examples with s_ = 

This completes what we will say by way of systematic analysis for relatively small N. 
In this subsection, we study another kind of example in which the factorization problem 
is easily solved explicitly. This is the case that all the double roots are in one factor of 
the M = 2 curve - for example, s_ = and therefore s+ = N — n. There is an analogous 
discussion for s + = 0. The condition < s + < N/2 shows that this is possible only when 
n < N < 2n. The factorization problem is solved by 



P N (x) +2A N = H 2 s+ (x)R N - 2s+ {x) 
P N (x) - 2A N = H 2 s+ {x)R N - 2s+ {x) - 4A 



N 



(3.63) 



with arbitrary coefficients in H 2 (x) and Rn-2 S+ (%)■ 



From ( |3.63j ), we easily find the semiclassical limit 



P N (x) -> H 2 (x)R N -2s + (x). 



(3.64) 



We see that U(N) is broken to U(2) s + x ^(l)^" 2 ^. 
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From ( |3.63| ) the matrix model curve is 

y 2 m = R N - 2s+ (x) [h? + (x)R N -2s + {x) -4A*] = [H s+ (x)R N - 2s+ (x)] 2 - 4A N R N - 2s+ (x). 

(3.65) 

If n + 1 < N, the degree of Rn-2 S+ (x) is N — 2s + = 2n — N < n — 1, and we 
immediately find the superpotential and /(#): 

W'(x) = H a+ (x)R N -2s+(x) 

(3.66) 

/(x)= -4A Ar J R Ar _ 2s+ (x) 

If n + 2 < N, the degree of f(x) is smaller than n — 1 and we conclude that S = 0. 

We have seen an example of this when we studied n = 2 with N = 4 above with 
(s+,s_) = (2,0). As a more interesting example, consider = 5 with (s + ,s_) = (2,0) 
and therefore n = 3. Following the general solution (|3.63|), we find 



P 5 (x) + 2A 5 = (x 2 + ax + b) 2 (x + c) 

(3.67) 

P 5 (x) - 2A N = (x 2 + ax + b) 2 (x + c) - 4A 5 . 

In the semiclassical limit, P§(x) — > (x 2 + ax + b) 2 (x + c), and therefore U(5i) is broken to 
U{2) x U{2) x U(l). 

From ( |3.65| ),( |3T66l ), we deduce 

Um= [(x 2 + ax + b)(x + c)] 2 -4A 5 (x + c) 
W'(x) = (x 2 + ax + b)(x + c) 

(3.68) 

f(x) = -4A 5 (x + c) 
S = 0. 

The theory has a quartic superpotential, and therefore there is no global symmetry. Yet 
we still find S = 0. This is a special case of a more general phenomenon we will see in the 
next section. 

3.7. Interpolation and duality 

As discussed in section 2.4, a given branch, Coulomb or confining, can have different 
classical limits. We have given explicit examples of the interpolation between classical 
limits U(Ni) x J7(AT 2 ) and U(Ni) x [/(A^). Since the A^ are certain periods of the one- 
form T, such an interpolation must involve a rearrangement of the compact one-cycles on 
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the Riemann surface. With some numerical work, we have determined how this happens 
in all the examples of section 3. Here we will simply illustrate the result qualitatively, 
without attempting any proofs. In the process, we will see that the two descriptions are 
related to each other through an electric-magnetic duality. 

For a general breaking pattern U(N) — > Yl7=i ^C^i)> there are 2n — 1 compact cycles, 
as sketched (for n = 3) in figures 1 and 2. There are n cycles Ai, . . . , A n that in a weak 
coupling limit surround a pair of roots, and n — 1 cycles Cj that connect neighboring pairs. 
The intersection numbers are Ai n Aj = dC\ Cj = 0, Ai n Cj = <$i,j+i — In particular, 
the sum i> = Ai+A2 + ...+ A n is a null vector. In fact, it is clear from figure 1 that this 
cycle is homologous to a cycle at infinity that winds once around the whole x-plane; this 
cycle will not be affected by the motion of the roots. So the null vector v will be invariant 
under the monodromy. If the action of the monodromies on Ai and Cj, i,j < n, is known, 
then the behavior of A n is determined using the fact that v is invariant. 

To illustrate, we will examine examples with n = 2. In this case, there are just three 
cycles, Ai, A 2 , and C = C\. If we know what happens to A\ and C, the behavior of A 2 is 
known. Given the intersection relations of A\ and C, the monodromies must act on them 
by an SL(2, Z) transformations plus possible addition of a multiple of the null vector v. 
(The null vector is uniquely determined by the fact that the Ni are no n- negative.) 

Now we will sketch the mechanism of interpolation we found to operate in all examples. 
These branches all have complex dimension two, but one complex parameter is associated 
with the freedom of shifting x in Pn(x). This parameter can be removed, leaving only 
a space with complex dimension one. We want to discuss how the zeroes of Pn(x) 2 — 
4A 2N move and the cycles evolve as one interpolates from a classical limit with unbroken 
U(N-l) x U(N 2 ) to a classical limit with U(Ni) x U(N 2 ). 

We start in the semiclassical limit with the group U(N\) x U(N 2 ). As discussed in 
section 2, the polynomial Pn(x) 2 — 4A 2N is given in this limit in terms of the Chebyshev 
polynomials of the two factors U(N 1 ) and U(N 2 ). Each factor leads to Ni — 1 double roots 
and two single roots. This is depicted in figure 4a for N± = N 2 = 2. 
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Figure 4: Interpolation between classical limits in the Coulomb branch of £7(4). The polynomial 
P%(x) — 4A 8 has four single roots (labeled 1,2,3 and 4) and two double roots; single and double 
roots are depicted by black and white circles respectively (a) £7 (2) x £7(2) semiclassical limit. 
Natural basis of compact cycles A\, A2 and C. (b) Semiclassical regime of £7(2) x £7(2) and 
transition to strong coupling of £7 (1) x £7(3). (c) Strong coupling of £7 (2) x £7(2) and transition to 
semiclassical regime of £7(1) x £7 (3). (d) £7(1) x £7 (3) semiclassical limit. Natural basis of compact 
cyles A\ , A2 and C 

As we vary the superpotential to go to a strongly coupled region, the two groups mix 
(figures 4b and 4c) and the zeroes then separate into two new groups (figure 4d). One 
group has two single roots and Ni double roots, and the other has two single roots and 
N 2 double roots. 

It is easy to see the way the cycles are transformed. Following the steps from figure 
4a to figure 4d, A\ evolves into — C and C evolves into A\. We can describe this by the 

51 



monodromy matrix 

In this example, with this choice of cycles, the monodromy is in SL(2, Z), with no addition 
of a null vector. In fact, the transformation matrix in ( |3.69|) is the S generator of SL(2, Z), 
revealing the electric-magnetic nature of the duality. Since A\ + A 2 = A\ + A2, it follows 
from the above that 

A 2 =A 1 + A 2 - C. (3.70) 
Since Ni = § A T and c = § C T, we have 

N x =c 

N 2=Nl +N 2 -c (3.71) 
c = -N v 

For example, if we start with the Coulomb vacuum N± = N 2 = 2, c = 1, then we get 
N± = 1, iV~2 = 3, c = —2. This is the case depicted in figure 4. For producing such 
an interpolation, we used Pa{x) in (|3.30 ). It is convenient to introduce a dimensionless 



parameter — . The figures, from 4a to 4d, show the roots of P^(x) 2 — 4A 8 as the parameter 
changes from to 00. Other examples from section 3 can be described similarly. 

Let us now interpret this in the M = 1 theory, where superpotentials of the form 
( 2.131) describe the low energy dynamics. To construct the Af = 1 effective superpotential 



for the gluino bilinears Si, we need the periods of the differential R(x) as well as those of 



T. We have as in g] and eqn. ( |2.8| ) 



Si = <b R, 

1 8T 7 ^ 
R. 



2ni dSi J B 

To evaluate these expressions, we need the transformation laws of the noncompact cycles 
E>i. These are uniquely determined modulo possible addition of a null vector by the 
intersections with the compact cycles. Since the monodromy must also preserve BiC\B 2 = 
0, it is fixed modulo addition of a common null vector to B\ and B 2 , and this is irrelevant 
in the sense that it is equivalent to a 27r shift in the theta angle. So we get 

B x = Ai + B 2 

(3.73) 

B 2 = B 2 . 
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Consequently, the transformation laws are 




(3.74) 



1 dT _ 1 dT 

Q I 

2m dS x 1 2ttz dS 2 

1 dT 1 dT 
2m dS 2 2%i dS z 



4. Equations Characterizing Components Of Moduli Space 

With the exception of section 3.6, we have discussed so far examples in which the 
rank n of the low energy gauge group is small and the number w = N — n of condensed 
monopoles is correspondingly large. Here we will concentrate on examples with relatively 
small w, and correspondingly a large rank of the low energy gauge group U(l) n . 

If U(N) is broken classically to YliU(Ni), then w = — 1)- This implies that 

if w < N/2, a typical range in the following discussion, then some of the iV; are 1. This 
means that all vacua are Coulomb vacua and confinement is not relevant in distinguishing 
different branches. Nevertheless, there are different branches of the moduli space of vacua, 
for a reason explained in section 2.4. The different branches correspond to decompositions 
w = s + + s_, where s + and s_ are, respectively, the number of double roots of the two 
polynomials P N (x) + 2A 2N and P N (x) - 2A 2N . 

Not having confinement as a useful order parameter and with only a very limited role 
for discrete symmetries, there appear to be no conventional order parameters that explain 
the existence of w + 1 branches labeled by s+. Instead, as sketched in the introduction, 
we will seek here to characterize each branch by describing holomorphic functions of chiral 
operators that have a vanishing expectation value on a given branch. The functions with 
this property will depend on s+ and s_. 

Certainly, there is no loss of essential generality in focusing on the case s + > s_. If 
s + = S- or s + = S- + 1, then the analysis below does not reveal any unusual equations 
obeyed on a given branch. For each w, precisely one of these possibilities is realized, and 
we know of no unusual equations obeyed on this branch. We will find such equations in 
all the other cases 



s+> s-+ 2, 



(4.1) 
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and these relations will depend on s_. 

To make explicit the double roots, we write 

, 2N = H 2 S 

.2N _ 5-2 



Pn + 2A = H s+ Rn-2 S+ 



(4.2) 



P N —2A — H s _R N _ 2s _. 

The curve of the M = 2 theory is 

y 2 = (P N + 2A 2N )(P N - 2A 2N ) = (H s _H s+ ) 2 R N - 2s+ R N _ 2s _. (4.3) 

However, in the matrix model we remove the quadratic factor (H S _H S+ ) 2 , setting y m = 
y/H s _H s+ , so 

y 2 m = R N - 2s+ R N -2 S _- (4.4) 



We write this 

(4.5) 



2 Rn-2s + [H 2 s+ R n _ 2s+ - AK 2N ) ( H s+ R N _ 2s+ \7, 4A 2 ^ 



tf 2 V / V h 2 + r n _ 2s+ 

As shown in [17], the order parameters £ r = — (l/327r 2 )Tr <& r W 2 are 



r H s+ R N . 2s+ i 4A 2 ^ 



t r =<b dxx r y m =(h dxx r — ^ — w 1 — . (4.6) 

ioo ioo if s _ y n s+ R N - 2s+ 

Here the contour is over a large circle near infinity in the complex x-plane (and we recall 
that a factor of l/2i:i is included in the definition of the symbol §). A standard argument 



shows that the factor y 1 — 4A 27V / H 2 + Rn- 2s+ can be replaced by 1 if 

r<s+ + S--2 = w- 2 (4.7) 

as then the terms obtained by expanding the square root vanish too rapidly at infinity to 
contribute. The number of values of r for which this inequality is obeyed is s+ + s_ — 1; 
given our assumption (|4.1[ ), this is at least 2s_ + 1. 
For r such that (|4.7|) is obeyed, we have simply 



I r r H s+ RN-2s + , . s 

t r =(b dxx r — ^ (4.8) 

J oo H s _ 

The rational function H s+ Rn- 2s+ /H s _ has poles only at the s_ zeroes of H a , so it has 
the form 

— S+ ^. JV 2s+ = — — h polynomial (4-9) 
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for some complex numbers Gtj, bj. Now we evaluate 



to give 



E 6 



7' 



(4.10) 



valid for r < w — 2. 

This equation expresses the + s_ — 1 quantities t r , < r < s + + s_ — 2, in terms 
of the 2s - quantities and bj. So the t r cannot be independent; they will obey at least 
s + — S- — 2 independent algebraic equations. Before working out these equations in general, 
let us consider the first few cases. If s_ = 0, we have simply t r = 0, < r < + s_ — 2. 
If S- = 1, we have t r = ba r for some a and 5, whence trfj = tkti whenever i + j = k + I 
(and i,j,k,l < w — 2). The first such relation (and the only example if s_ = 1 and 
s+ = S- + 2 = 3) is the equation ^2 — tf = 0. We can write this equation as the vanishing 
of the determinant of the 2x2 matrix 



M 9 . 



One way to prove that det M 2 = for s 

M 2 = 



to ti 
h t 2 



1 is to write 



(4.11) 



1 

a i 



1 

a 2 



bx 








0,1 
«2 



(4.12) 

where we have written a formula for M 2 that is valid if s_ = 2. But if we set s_ = 1, we 
should take b 2 = a 2 = (and then a± = a, b\ = b) and clearly det M 2 = since the second 
factor in M 2 has vanishing determinant. 

In general, for any n, we set M n to be the matrix whose i,j matrix element is ti+j- 2 . 
For example, the next case is 




M 3 



If S- = n, M n can be written as a product of three matrices, namely 



with 



A, 



(4.13) 





- A B A T 




( 1 


1 


1 \ 


ai 


a 2 










a\ 


a\ 





(4.14) 



.r 1 



.r 1 



,n-i 



(4.15) 
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and B n a diagonal matrix with eigenvalues &i, 62, • • • , &n- (Thus, the i, j matrix element of 
A n is a l ~ 1 . ) A n is a Vandermonde matrix, whose determinant is 



The determinant of M n is hence 



detM n = n( a *- a J-) 2 n> ( 4 - ir ) 

i<j k 

and in particular is generically nonzero for s- = n. 

det M n depends on the t r with r < 2n. Since our derivation has assumed that r < 
w — 2, the maximum value of n that we should consider is 

'1 



n 



w — 1 
2 



(4.18) 



where for any real x, [x] denotes the greatest integer less than or equal to x. In what 
follows we set n to this value. If s_ < n, the formula (|4.14| ) for M n remains valid, but we 
should set all but S- of the £>, (and a«) to zero. Hence the matrix B n , and consequently 
also M n , has rank s_. It follows that the rank s_ + 1 minors of M n (the determinants 
of (s- + 1) x (s_ + 1) matrices obtained by omitting some rows and columns from M n ) 
vanish. These are our relations. 

Each minor is a homogeneous polynomial in the tj of degree s_ + 1. For example, for 
S- = 0, the rank s_ + 1 minors are simply the tj themselves, and for s_ = 1, they are 
the quadratic functions ttfj — t^ti (for % + j = k + I) whose vanishing we found earlier for 
s- = 1. 

Since M n is an (n + 1) x (n + 1) matrix, it has minors of all ranks up to n + 1 (the 
rank n + 1 minor being simply det M n ). The relations we have found are hence non-trivial 
whenever n > s_, or equivalently whenver s+ > s_ + 2. For s+ = s_ + 2, we have n = s_, 
and the unique relation of this type is detM s _ = 0. 

For any given s_, the rank s_ minors do not generically vanish, since one of them (the 
determinant of the upper left s_ x s_ submatrix of M n ) is detM s __i, which is evaluated 
by setting n = S- in ( [4.17|) . So a branch with given s_, and s+ > s_ + 2, is characterized 
by saying that the rank s_ minors of M n are generically nonvanishing, while the rank 
S- + 1 minors vanish identically. We have accomplished our goal of finding holomorphic 
functions of chiral order parameters whose vanishing on some branches and not others 
distinguishes the different branches, at least modulo the exchange s+ <-» s_. 
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Appendix A. Generalization to Superpotentials of Arbitrary Degree. 

In section 2.2 we restricted the degree of the superpotential to be less than N + 1. 
The main reason was that only the first N coordinates u^s are independent. In this section 
we will generalize the discussion by allowing a superpotential of any degree. The idea is to 
regard all u^s as independent coordinates on a larger space subject to constraints. Clas- 
sically, a simple way of obtaining the constraints is by looking at the generating function 

of ut-S, 



1 _ y- Tr ¥ cl _N ^ lu t 

rf (f) , L-^i T T-' + l ' ' 

tXJ J- f->\ tkj tAj tAj 

cl 1=0 1=1 

Notice that this object is also written in terms of Pn(%) = det(x — $ c ^) which is a function 
only of u^s for i = 1, . . . , N as follows, 

Tr -J— = ^1 = 4-logP(x). (A.2) 
x-$ P N (x) dx 5 K ' K ' 



Integrating (|A.2| ) and (|A.1[ ) with respect to x, we get, 

/ oo 



i=l 



where the integration constant was fixed by matching the x N terms. 

The constraints we are after come from imposing (|A.3|) i.e, the Laurent series has to 
terminate at order x°. Since the coefficient of x~ l only depends on m with % < I and is 
linear in u\ , imposing that to vanish gives a system of equations for u\ with / > N + 1 in 
triangular form that can be solved in terms of the first N coordinates. 

In writing the effective superpotential (|2.26|) one has to use (Tr $ z ) instead of Tr & cl . 



In the first section we did not make any distinction since (Tr ¥} = Tr ¥ cl for / < 2N - 1. 
For I > 2N the relation receives instanton corrections. This implies that a modification to 
the constraints is needed. We now turn to this issue. 
Consider the quantum mechanical analog of (|A.1[), 



1 \ N ^ lUi „ N 

> = — + > -ttt (A.4) 

x - $ / x ^ x l+1 



where we have denoted j(Tr $ z ) by Ui. 
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The full quantum mechanical analog of ( A. 2 ) is 
1 \ _ P' N (x) 



Tr 



Integrating with respect to x both (|A.4|) and (|A.5|) we get 



^ log ^(x) + y/ptW-U?^ (A.5) 



P JV ( a ;) + Jp2 (x) _ 4A 2iv 



A' 



exp - J2 



Ui 



i=l 



X- 



(A.6) 



where the constant of integration was determined to be | by taking A to zero. 
Finally we can find the analog of ( |A.3|) by solving for Pn(x) in ( A.6 ), 



rw x N ^UA A 2N ^ 

P N {x)=x exp I - ^ — I + exp I ^ 



v i=l 



(A.7) 



Now we are ready to write the generalization of (|2.26|) for a superpotential of any degree. 

k 



W eff = Y,9rU r+1 + jV k - N {x) ^exp ^-^-ij +_exp ] ] dx 

N-n 



+ £ (^( / - +Bj 

i=1 \ J x Pi J 



Pn(x) 
(x -Pi) [ 



-dx 



(A.8) 

where Vfc-Ar(x) is a polynomial of degree k — N whose coefficients are to be thought of as 
Lagrange multipliers imposing constraints determining U\ for Z = iV+l,...,fc + lin terms 
of Ui with I = 1,...,N and A. 

The next step is to follow the computations in section 2. The derivative of W e ff with 
respect to U r +i is 



9r + 



+ 



Vk-N(x) ( N 

■x exp 



X 



r+1 



A 



27V 



i=l 



X 



N 



exp 



Qn-u-i(x) dP N (x) 
H N _ n (x) dU r+1 



(A.9) 



dx = . 



Using (|2.31| ) and the fact that in the relevant range of r we have U r = u r 

dP N (x) 



dUr+l { 



r + l<N 



However, inside the integral we can write -s 

the integral vanishes. 



r + 1 > N 

for any r since for r + 1 > N 



58 



In order to simplify ( |A.9| ) note that after imposing the Vk-N(x) constraints, 

xN ex p ( - E = \ ( Pn{ ^ + \/pUx)-^ 2N ) + o(x k+i ). 

It is easy to see that the 0(x k+1 ) terms do not contribute to the integrals. Therefore, we 
can write the integral in (|A.9| ) that contains Vk-N(x) as, 

2A 2K 



Pn(x) + V^') - 4A 2N 



dx 



Using these results we can multiply (|A.9|) by z r and sum over r from zero to k to get, 
Finally, using ( p. 24] ) we can replace y/P^(x) — 4A 2N by a/ F2n(x)HN- n (x) in the first 



integral and Pn(x) by ^^/F^Jx)HN-n(x) + 0(x~ N ) in the second. As in section 2.2, the 
0(x~ N ) does not contribute to the integral. Note that had we not included instanton 
corrections in the constraints, we would have obtained Pn(x) instead of a/ P^j(x) — 4A 2N 
in the first integral and the 0(x~ N ) could not have been dropped for k > 2N. 
The final result is thus 

W'{z) = j^F 2n {x) {V k - N {x)H N - n {x) +Q N - n - 1 {x)) T^T^y- ( A - 10 ) 

This agrees with the matrix model equation of motion 

W'(z)= I Vm ^ dx (A.ll) 
J x - z 

if the matrix model curve is 

vL( x ) = F 2n (x)Ql_ n (x). 

where Qk- n (x) = Vk-N{x)H^- n {x) + QN- n -i(x). Moreover, ( |A.10| ) implies that is 
known up to a polynomial fk-i(x) of degree k — 1, i.e. 

yl = F 2n (x)Ql_ n (x) = W' k \x) + 

providing the generalization of (|2.25| ) and the result ( |2.41| ) discussed in section 2.2. 
Let us consider some special cases: 

1. No massless monopoles n = N: Then Qk-N(x) = Vk-N{x). 

2. Degree of superpotential equal to N + 1, i.e. k = N: Qk- n (x) = VoHjsr- n (x) + 
Qn- u -i(x). In particular, for n = N, i.e. U(N) completely broken to U(l) , Qo is a 
constant and y^x) = F2n(x) = P%{x) — 4A? N . 
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Appendix B. Proof of the Generalized Konishi Anomaly from Strong Coupling 
Analysis. 

In this appendix, using the results from section 2.2 we show that the generalized 
Konishi anomaly equation 

follows from the effective superpotential ( |2.26[ ). For simplicity we will assume that the 
degree of W is less than 2N + 1 such that (Tr W'($)) = Tr W'($ c i). 

Instead of viewing u r as the coordinates on the M = 2 Coulomb branch, we can use 
$ c z and mod out by U(N). This is valid except at points where some of the eigenvalues of 
$ c z coincide. Then, instead of varying W e ff with respect to u r as in (|2.33| ), we vary with 
respect to 0j, the eigenvalues of <fr c i (recall that Pn(x) = Yli( x ~ to get 

(We have used the result P>i = from section 2.2.) From this we derive that with z outside 
the contour of integration 

^ W'(^) = fy £ / Pn(x) dx 

Z-<$>cl 1 J (Z - (f>l){x - (f)l)(x - Pi ) 

^ ' / (^-^(cc-pi) \X-(f)i Z-(piJ 



1=1 i=l 
N-n 



i=l 



V Li i f Tr l — Tr 1 —— ) dx 

^ J (z-x){x-pi) V x-$ c i Z-QdJ 



Using ( |2.27|) Pjv(Pi)Tr p .-$ ; = 0, and the fact that Pn(4>i) = the first term does not 
contribute. The second term can be simplified using ( |2.35|) 



(B.4) 

where we used manipulations as in ( |2.36| ). It is important to keep in mind that we can 
not replace Pn(x) by a/ F2n( x )H n _ n (x) as we did in ( |2.37|) . The reason in that the terms 
of order 0(x~ N ) can not be dropped in the integral. To see this, note that \z\ > \x\ 
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for any x inside the contour of integration. Therefore, l/(z — x) can be expanded as 
a/z)J2Z (x/z) 1 . 

A convenient way to deal with this is to write, 



P N {x)Qk-n{x) 

out H N - n (x)(z - x) 



dx 



PN{x)Qk-n{x) 
H N - n (x)(z - X) 



dx 



PN(x)Qk-n(x) 

c H N - n (x)(z - x) 



dx 



where C z is a small contour around z. u out" and "m" refer to the point x = z being 
outside or inside the contour of integration. The first integral on the rhs gives W'(z) as 
( p.36| ) indicates. The second one can be evaluated at the pole. 
Using (|2.24|) to write, 



H 



N-n 



y/P*(z)-4A™ 



and y m (z) from (|2.39|) , we get ( |B.4|) to be 



W'(z)-y m (z) 



Pn(z) 



y/P*(z)-4A™ 



;b.5) 



The left hand side becomes a polynomial in z when combined with the term proportional 
to W'{z) on the right hand side. The right hand side can also be simplified by using that 



Tr 



1 = P'nW 
2-* ci Pn(z) ' 



(|B.5|) is then, 



In the right hand side we use the relation in the chiral ring of M = 2 



(B.6) 



Tr 



1 



a/ P'n(z) — 4A 2N \ ' z-<P 
For the left hand side we use the fact that 

W'($ cl )-W'(z 



Tr 



z - $d 

contains only Tr $ cl with I < k which for k < 2N is equal to (Tr $ ), 



Z — $ 
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Finally, we can write ( |B.6| ) as follows, 



Tr > = ( Tr -±- > (W(z) - y m (z)) (B.7) 



z — <& / \ z - <E> 

By definition, the resolvent of the matrix model -R(z) satisfies 2i?(z) = W'(z) — y m (z) 
Therefore ( |B.7| ) becomes the generalized Konishi anomaly equation, 

Tr-A_i)=2^)(Tr 



z — / \ z - $ 



Appendix C. Chebyshev polynomials 

In this appendix we list the first six Chebyshev polynomials of first and second kind. 
These are used in the examples discussed in section 3. 

The definition given in section 2.3 of Chebyshev polynomials of the first and second 
kind of degree t and t — 1 respectively is the following. Set x = cos 9 and let 

s f r\\ , , ! \ 1 d% , \ sin(t6>) 

%{x) = cos ^ U t -i(x) = --^{x) = —V- C- 1 

t ax smv 

From this it is simple to compute the first of then, 

t %(x) U t (x) 

1 x 2x 

2 2x 2 - 1 Ax 2 - 1 

3 4x 3 - 3x 8x s - 4x 

4 8x 4 - 8x 2 + 1 16x 4 - 12x 2 - 1 

5 16x 5 - 20x 3 - 5x 32x 5 - 32x 3 + 6x 

6 32x 6 - 48x 4 + 18x 2 - 1 64x 6 - 80x 4 + 24x 2 - 1 



Appendix D. A Magnetic Index 

Here we want to explain an interesting result that was obtained in an unsuccessful 
attempt to find a new order parameter for these models. The attempt was unsuccessful 
because the object v that we defined turned out to be computable in terms of objects that 
are already known. 

We will examine closely the physical spectrum of particles in a theory in which clas- 
sically U(N) is broken to a low energy subgroup such as U(Ni) x U(N2). (Eventually 
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we will generalize to an arbitrary breaking pattern U(N) — > ]X £7(iVj).) Allowing for the 
strong quantum dynamics, the true low energy gauge group is £7(1) x U(l). One £7(1) is 
the center of U (N) and decouples from the theory. We will disregard it and focus on the 
£7(1) that is not decoupled, whose generator we call Q. 

Physical particles in the theory are in general electrically and magnetically charged 
with respect to this £7(1). The product of the smallest nonzero electric charge q (of an 
unconfined magnetically neutral particle) times the smallest physical magnetic charge m 
is, according to Dirac, an integer v. (We measure m as a multiple of the Dirac quantum 
2nhc.) We will call this integer the magnetic index, and we claim that always 

vt = N, (D.l) 

where t is the confinement index defined in section 2.1. 

First let us verify this in examples. We consider first the case of £7(4) broken to 
£7(3) x £7(1). We embed £7(3) as the upper left 3x3 block in £7(4). For the monopole, 
we consider an 't Hooft-Polyakov monopole embedded in £7(4) in such a way that (in a 
unitary gauge) the magnetic charge is 

( \ 

° ! , (D.2) 

V -i) 

where in other words the "1" acts in £7(3) and the "—1" in £7(1). To find the smallest 
electric charge that is not confined, we start with the massive VF-bosons, which are a £7(3) 
triplet embedded as follows in £7 (4) : 

/0 W 1 

w 2 
w 3 

V 

The components W 1 , W 2 , and W 3 have, respectively, electric charges 1,1, and 2, with 
respect to the electric charge generator given in ( p.2| ).lll The SU(3) singlet combination 
€ijkW l W : 'W k is unconfined. Its electric charge isg = l + l + 2 = 4. Since the minimum 

12 In evaluating the electric charge of an SU(3) singlet, we can simply use the charge matrix 
that appears in ( |D.2| ); subtracting an SU(3) generator to get a multiple of Q would not change 
the result. 
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magnetic charge m in this theory is 1 (in Dirac units), the product v = qm equals 4. The 
breaking of £7(4) to £7(3) x £7(1) gives a Coulomb branch, so t = 1, and vt = 4, as claimed 
in O- 

Now we consider the breaking of £7(4) to £7(2) x £7(2). We embed the two factors of 
£7(2) as the upper left and lower right blocks in £7(4), and the monopole charge as 



1 

-1 

V o/ 



(D.4) 



We focus on the upper right 2x2 block of W bosons, transforming as (2, 2) of £7(2) x £7(2): 



W 11 W 12 



W 21 W 22 



(D.5) 



The components (W , W 21 , W 12 , W 22 ) have (with respect to the charge generator in 
( p.4| )) respectively electric charges (1,2,0,1), and the basic SU(2) x ££7(2) singlet 
6ij€i'j'W u = W 1X W 22 — W 12 W 21 hence has charge q = 2. In the confining branch, 
with b = in the notation of section 2.1, the confining index is t = 2, the minimum 
magnetic charge is again m = 1, and v = qm = 2. So again vt = 4, as claimed. 

What, however, happens in the branch with classical breaking to £7(2) x £7(2) and 
with 6 = 1? In this case, there is no confinement, and t = 1. What about vl Shifting 
the value of b changes us from confinement to oblique confinement in one of the SU(2) 's, 
but this does not affect which combinations of W bosons are confined. So the minimum 
electric charge, in the same units, is still q = 2. However, in section 2.1, we described a 
process of magnetic screening that occurs in this model when an external Wilson loop in 
the fundamental representation of SU (4) is considered; it is screened by the nucleation 
of an 't Hooft-Polyakov monopole with unit charge, and this is the mechanism that leads 
to t = 1. The fact that the 't Hooft-Polyakov monopole screens an external Wilson line 
means that, in the absence of the Wilson line, it would have infinite energy and would be 
confined (we explain the mechanism in more detail below). So the minimum unconfined 
magnetic charge when 6=1 actually has m = 2. Since q is still 2, this gives v = 4, which, 
with t = 1, is again consistent with vt = 4. 

At least in hindsight, it should not come as too much of a surprise that we had to 
consider here the confinement of certain magnetic monopoles. We have taken account of 
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electric confinement to evaluate q, so we should expect to consider magnetic confinement, 
that is confinement of some magnetic monopoles, to evaluate m. 

More directly, the confinement for b = 1 of the minimum charge monopole arises as 
follows. The monopole has a collective coordinate whose quantization leads to the Julia- 
Zee dyon. There is only a single collective coordinate, and quantizing it gives the monopole 
equal color charges in the two S'C/(2)'s£l In SU(2) gauge theory, ordinary confinement 
means that a monopole is unconfined if its color electric charge is an integer multiple of the 
W boson color electric charge, and oblique confinement means that the color electric charge 
of an unconfined monopole is a half-integral multiple of the W boson color electric charge. 
Setting 6 = 1 means that one SU(2) has ordinary confinement and the other has oblique 
confinement. Given this, the states obtained by quantizing the collective coordinate are 
all confined, as we have claimed. 

Proof Of The Relation 

Now that we have given an idea of what the relation vt = N means semiclassically, 
we proceed to a mathematical proof of this relation. In the process, we will also define the 
index v for arbitrary rank of the low energy group; above, we assumed rank 1. 

We start with the M = 2 Coulomb branch. The low energy gauge group (omitting 
the center of U(N)) is U(1) N ~ 1 . The massive gauge bosons and monopoles have charges 
that fill out a rank 2(N — 1) lattice that we will call L. On L there is an integer-valued 
antisymmetric inner product that we will call ui. For example, if iV — 1 = 1, there is a 
single U(l) and a vector in L is specified by giving electric and magnetic charges (e,g). 
Given two vectors I = (e,g) and I' = (e',g'), the inner product is then u(l, I') = eg' — e'g. 
In general, for any number of £7(1) 's, ui is defined by a similar formula, summing over all 
of the C/(l)'s. The Pfaffian of uj is in absolute value 

|PfaffH| = N. (D.6) 

This relation reflects the fact that if we incorporate electric charges in the fundamental 
representation of SU(N), then uj would become unimodular. 



13 Global SU(2) charge is not defined in the field of a magnetic monopole [22 1. What we here 
call the color electric charge (for a given SU(2)) is the quantity that is well-defined. It originates 
from the projection of the charge generator in ( p.4|) into the SU{2) of interest; this charge is an 
integer for SU (2) generators, and would be a half-integer for external quarks. 
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Now, after perturbing the theory by a superpotential (or in any other way, for that 
matter), some charges may condense. The condensed charges generate a sublattice M of 
L. The condensed objects are always "mutually local," in the sense that for m, m! G M, 

u(m, m) = 0. (D.7) 

A particle is confined if its charge vector / is such that u(l, m) 7^ for some m G M. 
So the charge vectors of unconfined charges lie in M , the sublattice of L consisting of I 
such that u(l,m) = for all m G M. Physically, the charge of an unconfined particle can 
only be measured modulo M, since particles with quantum numbers in M have condensed. 
So the measureable charges in the low energy theory take values in L' = M /M. The 
antisymmetric product uj induces an analogous antisymmetric product u' for the physical, 
unconfined charges in the low energy theory. The general definition of v is 

v= |Pfaff(a/)|. (D.8) 

v determines how many extra types of charge could be added to the low energy theory 
(keeping fixed its gauge group) without any inconsistency. 

V is the lattice of charges that are not confined, and hence the quotient L" = L/L' 
might be regarded as the lattice of charges that are confined. There is a well-defined inner 
product uj" : M x L" — > Z; for any m G M, I" G L" , we simply lift I" to an element I G L 
in an arbitrary fashion, and define uj"(m,l") = uj(m,l). The choice of I is well-defined 
modulo an element of M -1 ; adding such an element to / does not change u(m, I). We can 
extend uj" to an antisymmetric form on M © L" by defining uj"(l" , m) = —oj(m, I"). The 
confinement index t is t = Pfaff(o/'). In fact, the confinement index measures how many 
types of charge are conceivable (they have integer inner products with vectors in M) but 
are absent in the theory and can be represented instead by external Wilson loops. 

Now it is a fact of linear algebra that in this situation 

Pfaff(cj) = Pfaff(u/) ■ Pfaff(u;"). (D.9) 

Putting all of this together, we get the promised relation iV = vt. 

If we are willing to make some unnatural choices, we can define the lattices and explain 
the key fact (|D.9| ) in a way that might be easier to follow. There is not any natural way 
to embed V as a sublattice of L, because vectors in V can only be interpreted as elements 
of L modulo M. However, this indeterminacy is unimportant for us; we simply embed V 
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as a sublattice L\ C L by picking a basis of V and lifting each basis vector to a vector in 
L, making a choice that (for each basis vector) is unique up to adding an element of M. 
If we do this, then u' can be identified with u±, the restriction of ui to L\. 

Having identified L\ as a sublattice of L, we can without further arbitrary choices 
define another sublattice L2 = L^; that is, L2 consists of elements 1% G L such that 
w(h,l2) = for all l\ G L\. With respect to the decomposition L = L\ © L2, uj is 
block-diagonal, 




(D.10) 



From this, ( p.9|) is clear. 
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